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IMPLEMENTATION  OF  THE  MARJORIZATION  PROGRAM 
FOR  COMPUTING  OPTIMUM  COMBINED  KINEMATIC 
VARIABLES  ON  THE  HP9000/835 


Introduction 

The  purpose  of  this  report  is  to  document  the  implementation  of  an  improved  and  expanded  version 
of  the  “Marjorization”  computer  program  on  the  HP9000/835,  and  to  present  detailed  instructions  on  the 
execution  of  this  program. 

An  underlying  assumption  of  the  Marjorization  process  is  that  photo-derived  displacements  output  by 
the  Naval  Biodynamics  Laboratory  s  data  analysis  program,  EZFLOW,  are  the  best  source  for 
displacement  variables,  and  EZFLOW  sensor-derived  accelerations  are  the  best  source  for  acceleration 
variables.  Another  assumption  is  that  small  corrections  to  various  accelerometer  constants  can  be 
determined,  which  produce  very  little  change  in  the  kinematic  variables  at  the  acceleration  level,  but 
which  can  produce  significant  changes  in  the  computed  displacements  and  can  bring  the  corrected  sensor- 
derived  displacements  into  agreement  with  the  photo-derived  displacements.  Marjorized  variables  are 
defined  as  the  corrected  accelerometer-derived  kinematic  variables  obtained  by  applying  the  computed 
corrections  to  the  appropriate  parameters  in  the  EZFLOW  sensor  program.  These  variables  represent 
the  optimum  combination  of  the  EZFLOW  photo-  and  sensor-derived  variables,  and  they  provide  a 
consistent  set  of  kinematic  variables  (from  the  acceleration  level  to  the  displacement  level)  suitable  for 
mathematical  modeling. 

As  discussed  in  detail  in  previous  publications,’’^  the  purpose  of  the  Marjorization  program  is  to 
determine  the  most  probable  corrections  to  accelerometer  direction  cosines  and  sensitivities  and  to  initial 
Euler  angles,  which  would  bring  the  accelerometer-derived  linear  displacements  and  angular  orientations 
into  agreement  with  the  cprresponding  photo-derived  variables.  The  Marjorization  algorithm  constrains 
the  solution  for  the  corrections  in  order  to  deter  error-component  solutions  that  excessively  exceed  the 
a  priori  values  for  these  components,  since  large  corrections  could  significantly  change  the  computed 
acceleration  components.  Appendix  A  includes  a  summary  of  the  Marjorization  process. 


Description  of  the  Marjorization  Program 

The  new  Marjorization  program  consists  of  three  primary  modules  that  may  be  executed  using  a  single 
UNIX  shell  script  xqtmarj  (Appendix  Bl).  A  flowchart  of  the  program,  using  the  input  and  output  data 
files  required  for  sled  run  LX3959,  is  provided  in  Figure  C-1.  The  data  files  for  the  same  run  are 
summarized  in  Table  D-1,  which  also  cites  the  tables  that  contain  descriptions  of  the  required  data  for 
the  various  modules  and  the  figures  that  contain  sample  data  files.  Appendix  B2  presents  sample  UNIX 
shell  scripts  that  may  be  used  to  execute  the  individual  Marjorization  modules  and  the  auxiliary  modules. 

The  first  module  computes  perturbed  displacement  variables  using  perturbation  values  read  from  a  data 
file.  The  second  module  uses  these  perturbed  variables  to  compute  partial  derivatives  for  each  of  the 
sensor  parameters  and  computes  a  constrained  least-squares  solution  for  the  corrections  to  the  sensor 
parameters  that  minimize  the  differences  between  the  photo-derived  and  sensor-derived  displacements 
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within  the  context  of  the  previously  indicated  constraints.  The  third  module  reads  these  corrections, 
applies  them  to  the  sensor  parameters,  and  computes  corrected  kinematic  variables,  which  are  defined 
as  Marjorized  variables.  There  is  also  an  auxiliary  module,  which  rewrites  the  reduced  EZFLOW  photo- 
derived  variables  in  the  format  required  by  the  second  Marjorization  module,  and  a  plotting  module, 
which  employs  the  standard  routines  used  at  the  Naval  Biodynamics  Laboratory  to  generate  EZFLOW 
plots. 

The  primary  modules  incorporate  several  significant  improvements  over  the  previous  version  of  the 
Marjorization  program,  which  was  implemented  on  the  UNIVAC  in  the  1980s.  Most  of  the  important 
changes  are  summarized  in  Table  1 . 


Table  1.  Major  Changes  in  Marjorization  Program. 

Language 

m 

EZFLOW 

Version 

Iterative 

Solutions 

Max  U  of  Data 
Pts  Used  in 
Sol'n 

Automated 

Solution 

Old 

Fortran  IV 

6 

Pre-1986 

No 

60 

No 

New 

Fortran  77 

9 

Current 

Yes 

150 

Yes 

The  code  is  now  written  in  FORTRAN  77  and  includes  all  of  the  changes  and  improvements  that  have 
been  implemented  in  EZFLOW  up  to  the  present  time.  In  addition,  the  Marjorization  program  now 
computes  corrections  to  all  available  accelerometers;  that  is,  corrections  are  computed  for  all  nine  head 
accelerometers,  while  in  the  past  calculations  were  limited  to  six  accelerometers,  since  the  original 
Marjorization  code  was  developed  when  only  six  were  used  on  the  head.  In  addition,  iterative  solutions 
may  now  be  obtained,  with  the  iteration  implemented  automatically  in  a  shell  script  (Appendix  Bl)  for 
a  predetermined  number  of  iterations.  The  corrections  computed  for  each  iteration  are  used  to  update 
the  reference  values  used  for  the  next  iteration,  and  each  successive  solution  is  linearized  around  the 
solution  from  the  previous  iteration. 

The  capability  of  inputting  factors  to  be  applied  to  the  a  priori  variances  for  the  individual  photo- 
derived  linear  and  angular  displacement  variables  has  also  been  implemented  in  the  new  Marjorization 
program.  The  manipulation  of  these  factors  was  found  to  be  the  most  effective  way  of  improving  the 
agreement  between  the  Marjorized  displacement  variables  and  the  photo-derived  displacements,  when  the 
Marjorized  variables  obtained  using  the  baseline  values  for  these  factors  produced  less  than  the  desired 
degree  of  agreement.  To  facilitate  the  use  of  the  Marjorization  program  with  Optotrak™  solutions,  the 
maximum  number  of  data  points  to  be  used  in  the  solution  has  also  been  increased  from  60  to  150. 


Results  and  Validation 

The  new  Marjorization  program  was  tested  on  thirteen  sled  runs  (Table  D-2).  Sample  plots  comparing 
the  Marjorized  variables  for  sled  run  LX3959  with  the  photo-derived  and  accelerometer-derived  variables, 
which  will  be  analyzed  later  in  this  paper,  are  included  in  Figure  C-2  through  Figure  C-13.  The  new 
program  was  exercised  extensively  on  some  of  the  runs,  and  based  on  these  results,  some  refinements 
were  made  to  the  previously  established  nominal  values  for  the  variable  input  data  values.  The  weighting 
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scheme  for  selecting  data  points  to  be  included  in  the  solution  was  modified  to  take  advantage  of  the  new 
capability  of  using  up  to  150  data  points  in  the  solution.  Taking  advantage  of  the  new  capability  to  apply 
factors  to  the  variances  for  the  linear  and  angular  displacement  variables,  baseline  factors  for  the  head 
and  T-1  were  determined.  In  conjunction  with  this,  refinements  were  also  made  to  the  factors  previously 
applied  to  the  a  priori  variances  for  the  accelerometer  direction  cosines  and  sensitivities,  as  well  as  to  the 
initial  Euler  angles. 

The  matrix-inversion  technique  used  to  obtain  the  least  squares  solution  for  the  error  components  was 
also  verified  by  implementing  a  method  for  obtaining  solutions  to  machine  accuracy  obtained  from  the 
NSWCLibrary  of  Mathematics  Subroutines,^  ''Solution  of  Equations  with  Iterative  Improvement.” 

A  comparison  of  the  solutions  based  on  the  two  methods  revealed  negligible  differences. 

Results  obtained  from  iterative  solutions  were  analyzed  and  the  following  observations  were  made. 
As  previously  reported,  it  was  found  that  a  single  iteration  generally  produced  satisfactory  agreement 
between  the  Marjorized  variables  and  the  photo-derived  displacements  without  excessively  large 
corrections  to  the  accelerometer  parameters.  However,  although  most  runs  were  handled  well  with  a 
single  iteration,  in  some  cases  enough  nonlinearity  exists  to  warrant  iterative  solutions.* 

Although  the  bulk  of  the  possible  corrections  to  the  sensor-derived  linear  and  angular  displacements 
are  obtained  by  applying  first-iteration  corrections  to  the  sensor  parameters,  it  was  found  that  the 
corrections  computed  for  two  iterations  generally  produced  even  better  overall  agreement  between  the 
photo-derived  and  accelerometer-derived  displacements,  without  resulting  in  excessively  large  corrections 
to  the  accelerometer  parameters.  It  was,  however,  generally  found  that  very  little  improvement  can  be 
obtained  with  more  than  two  iterations,  and  this  usually  comes  at  the  expense  of  larger  cumulative 
computed  corrections,  which  can  become  unacceptably  large  compared  to  the  a  priori  values  for  the 
errors  in  the  accelerometer  parameters.  In  addition,  further  iterations  generally  improve  the  off-axis 
displacement  components,  in  some  cases  at  the  expense  of  small  degradation  in  the  midsagittal-plane 
displacement  components,  particularly  the  Z  displacement. 

Recommended  Baseline  Data  Values.  Based  on  an  analysis  of  the  results  of  a  large  number  of  computer 
runs  for  the  thirteen  sled  acceleration  tests  used  in  this  study,  the  following  is  recommended.  Production 
solutions  should  be  obtained  using  two  iterations  and  applying  the  factor  sets  shown  in  Table  2  to  the 
denominators  of  the  photo-derived  linear  displacements  and  Euler  angles.  In  addition,  a  factor  of  1.0 
should  be  applied  to  the  variances  for  the  accelerometer  direction  cosines  and  sensitivities,  as  well 

as  to  the  initial  Euler  angles,  as  compared  to  a  factor  of  3.0  that  was  previously  used.  These  values, 
together  with  the  revised  weighting  scheme  for  the  selection  of  data  points  to  be  used  in  the  solution,  are 
included  in  the  sample  data  files  in  Appendix  B3. 


Table  2.  Factor  Sets  to  be  Applied  to  the  Denominators  of  the  Photo-Derived  Linear 

Displacements  and  Euler  Angles  to  Obtain  Production  Solutions. 

Linear 

Angular 

Head 

1 

10 

T1 

3 

3 

These  recommendations,  together  with  the  previously  established  nominal  values  for  Marjorization 
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input  data/  should  be  used  for  production  Marjorization  runs.  Nominal  values  for  all  user-required  input 
data  to  the  Marjorization  modules  are  shown  in  sample  data  files  provided  in  Appendices  B3  through  B8; 
Tables  D-3  through  D-8  contain  descriptions  of  this  data.  Table  D-9  defines  the  variables  included  in 
plot  specification  module /i/iwja/y.  Tables  D-10  and  D-1 1  contain  descriptions  of  program-generated  data 
that  is  input  to  follow-on  Marjorization  modules.  Sample  output  from  the  Photo-Reformatting  Module 
is  show  in  Appendix  B9.  In  addition,  Appendix  B1  provides  a  sample  shell  script  which  automatically 
cycles  through  the  desired  number  of  iterations  (in  this  case,  two),  deleting  data  files  as  necessary  in 
accordance  with  program  requirements.  All  of  the  execution  scripts  and  data  files  are  stored  on  the 
HP9000  in  the  dir^tory  /7933/users/margie/MARJORIZATION,  with  the  exception  of  the  plot-generation 
scripts,  which  are  contained  in  the  directory  /7933/users/margie/MARJORIZATION/PLOTS. 

Additional  iterations  and/or  modification  of  the  nominal  input  data  values  could  be  appropriate  if 
agreement  between  the  Marjorized  displacement  variables  and  the  photo-derived  variables  is  less  than 
desired,  or  if  the  computed  corrections  are  larger  than  desired  for  a  particular  sled  run.  However,  past 
experience  has  shown  that  less  than  satisfactory  results  from  the  Marjorization  process  often  is  an 
indication  of  some  problem  with  the  run  data,  which  should  be  investigated  and  corrected  if  possible. 
In  the  event  that  it  is  necessary  to  improve  the  Marjorization  results,  manipulation  of  factors  applied  to 
the  denominators  of  the  a  priori  variances  for  photo-derived  displacements  was  found  to  be  effective. 
These  factors  are  entered  in  the  input  file  fllein  (Appendices  B3  and  B6). 

Most  of  the  sled  runs  included  in  this  study  appear  to  have  a  noticeable  timing  discrepancy  between 
the  photo-derived  and  sensor-derived  X  displacements  for  both  the  head  and  neck.  Since  this  timing 
discrepancy  is  not  generally  seen  in  the  other  displacement  variables,  this  problem  appears  to  be  the  result 
of  a  discrepancy  between  the  computed  times  for  the  sled  variables  and  those  for  the  head  and  neck 
variables,  since  sled  displacement  is  subtracted  from  the  head  and  neck  X  displacements  in  the  laboratory 
coordinate  system  to  obtain  the  X  displacement  in  the  sled  coordinate  system,  which  is  the  plotted 
variable. 

For  example,  in  Figure  C-14  and  Figure  C-15,  in  which  the  top  plot  displays  the  head  X  displacement 
and  the  bottom  plot  the  T1  X  displacement,  it  may  be  seen  that  there  is  a  significant  discrepancy  in  the 
timing  of  the  sensor-derived  and  photo-derived  X  displacements  for  both  the  head  and  neck  for  runs 
LX3965  and  LX3904  respectively.  In  these  plots  the  last  character  of  the  Marjorized  variable  name  is 
C,  the  last  character  of  the  photo-derived  variable  name  is  P,  and  the  last  character  of  the  sensor-derived 
variable  name  is  S.  Since  timing  discrepancies  of  this  type  cannot  be  corrected  to  any  significant  degree 
with  normal  Marjorization  techniques,  the  Marjorization  program  was  modified  to  permit  shifting  the  sled 
acceleration  trace  in  an  effort  to  correct  this  discrepancy.  A  shift  of  only  1  msec  was  found  to 
significantly  improve  the  agreement  between  the  X  Marjorized  and  photo-derived  displacement  variables, 
as  can  be  seen  in  Figure  C-16  through  Figure  C-19,  in  which  the  top  plot  reflects  the  unshifted  sled 
acceleration  trace  and  the  bottom  one  a  1-msec  shift.  Note  that  although  there  is  still  a  significant 
discrepancy  in  the  timing  of  the  curves,  this  small  timing  shift  resulted  in  significantly  improved 
agreement  between  the  Marjorized  and  the  photo-derived  X  displacements  for  both  the  head  and  neck. 

Criteria  for  Evaluating  the  Marjorization  Results.  Objective  criteria  were  established  for  evaluating 
the  Marjorization  results  on  the  basis  of  analysis  of  a  large  number  of  runs.  These  criteria  are  displayed 
at  the  bottom  of  Tables  D-2,  D-12,  and  D-13.  A  Marjorized  linear  displacement  variable  was  Judged 
as  excellent  if  there  was  less  than  a  1.25-cm  difference  between  the  photo-derived  and  Marjorized  variable 
in  the  time  window  of  interest  (between  the  variable  peak  and  250  msec).  An  angular  displacement 
variable  was  Judged  as  excellent  if  there  was  less  than  a  three-degree  difference  between  the  photo-derived 
and  Marjorized  variable  in  the  time  window  of  interest.  Because  acceptable  Marjorization  results  also 
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require  that  the  corrections  to  the  accelerometer  parameters  not  be  excessively  large  for  runs  without 
known  problems,  evaluation  criteria  were  also  established  for  the  corrections  to  the  accelerometer 
orientations  and  sensitivities,  as  well  as  corrections  to  the  initial  Euler  angles. 


An  Analysis  of  Results  for  Run  LX3959 

Plotted  results  for  run  LX3959  comparing  accelerometer-derived  displacement  variables  (last  character 
of  the  variable  name  S),  photo-derived  variables  (last  character  P),  and  Marjorized  variables  (last 
character  Q  are  presented  in  Figure  C-2  through  Figure  C-13.  Although  these  plots  were  generated 
before  the  program  modification  that  increased  the  number  of  data  points  used  in  the  solution  from  60 
to  150,  these  results  do  not  differ  significantly  from  those  obtained  with  the  new  baseline  input  data 
values.  The  following  analysis  is  based  on  the  data  in  Appendix  B6.  However,  a  similar  analysis  would 
apply  using  the  modified  input  data  in  Appendix  B3.  Computed  corrections  that  were  applied  to  the 
sensor  parameters  to  generate  the  Marjorized  variables  are  shown  in  Table  D-14. 

It  may  be  seen  from  Figure  C-2  through  Figure  C-7,  for  which  the  plotted  variables  are  defined  in 
Table  D-9,  that  while  the  Marjorized  curves  for  the  head  displacement  variables  closely  approximate  the 
photo-derived  curves,  the  sensor-derived  displacement  variables  drift  significantly  from  the  photo 
solution,  with  the  exception  of  the  second  Euler  angle  in  Figure  C-6,  in  which  the  sensor-derived  curve 
reasonably  approximates  the  photo-derived  curve.  In  addition,  although  the  Marjorized  neck  displacement 
variables  shown  in  Figure  C-7  through  Figure  C-13  do  not  agree  as  well  with  the  photo-derived  variable 
as  do  the  head  variables,  the  magnitude  of  the  differences  in  the  time  window  of  interest  are  quite  small, 
while  the  sensor-derived  neck  displacement  variables  all  drift  away  from  the  photo  solution. 
Furthermore,  as  may  be  observed  in  Figure  C-8,  there  is  an  obvious  timing  discrepancy  between  the 
photo-derived  and  sensor-derived  X  displacement  variables,  and  a  small  shift  in  the  time  of  the  sled 
acceleration  would  be  expected  to  improve  the  neck  X  results,  as  previously  discussed. 

We  will  now  apply  a  more  objective  analysis  to  Figure  C-2  through  Figure  C-13  using  the  evaluation 
criteria  shown  at  the  bottom  of  Table  D-2.  As  may  be  observed  from  the  plots,  with  the  exception  of 
the  neck  X-displacement,  all  head  and  neck  Marjorized  variables  differed  from  the  photo-derived  variables 
by  less  than  1.25  cm  during  the  time  window  of  interest  (between  the  curve  peaks  and  250  msec)  and 
would  therefore  be  categorized  as  excellent.  The  neck  X  displacement  exceeds  1.25  cm  by  only  a  small 
amount,  and  therefore  is  classified  as  G-t-,  i.e.,  good  plus.  In  addition,  all  of  the  Marjorized  head  and 
neck  Euler  angles  differed  from  the  corresponding  photo-derived  Euler  angles  by  less  than  3°  and  were 
therefore  classified  as  E  (excellent).  The  magnitudes  of  all  of  the  computed  corrections  to  the  head  sensor 
parameters  (Table  D-14)  were  classified  as  excellent,  although  the  correction  to  the  sensitivity  for 
accelerometer  5  was  —3.5%,  which  is  the  lower  limit  for  the  excellent  category.  Hence  the  overall  rating 
for  the  “Magnitude  of  Corrections”  is  E— .  Because  the  computed  correction  to  the  orientation  of  neck 
accelerometer  2  was  4.75  degrees,  which  is  classified  as  good  according  to  the  evaluation  criteria,  the 
overall  “Magnitude  of  Corrections”  rating  for  the  neck  is  G.  Although  the  neck  X  displacement  was 
classified  as  G-I-,  a  small  shift  in  the  sled  accelerometer  trace  would  be  expected  to  bring  the  Marjorized 
X  displacement  into  better  agreement  with  the  photo-derived  variable  and  this  variable  could  then 
presumably  be  categorized  as  excellent. 

Hence  it  may  be  seen  that  the  small  corrections  to  the  accelerometer  parameters  and  initial  Euler  angles 
shown  in  Table  D-14  resulted  in  significant  changes  in  the  corrected  accelerometer-derived  displacement 
variables,  that  is,  the  Marjorized  variables.  Although  no  effort  was  made  to  further  improve  the 
agreement  between  the  Marjorized  and  photo-derived  neck  X  displacement  by  manipulating  input  data 
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values,  it  is  believed  that  any  substantial  improvement  in  this  variable  would  require  shifting  the  sled 
acceleration  trace. 


An  Analysis  of  the  Marjorized  Results  for  Other  Sled  Runs 

Table  D-2  presents  an  evaluation  of  the  Marjorized  displacement  variables  based  on  new  baseline  input 
data  values  for  all  thirteen  sled  acceleration  tests  included  in  this  study.  Both  the  head  and  the  neck 
displacement  variables  were  categorized  according  to  the  criteria  displayed  at  the  bottom  of  the  table. 
It  may  be  seen  that  for  9  of  the  13  runs,  all  head  Marjorized  displacement  variables  were  judged  to  be 
excellent,  and  with  the  exception  of  runs  LX3959  and  LX3970,  the  magnitude  of  the  corrections  to  the 
head  accelerometers  for  each  of  these  runs  was  also  classified  as  excellent.  For  the  neck,  no  linear  or 
angular  displacement  variable  was  categorized  as  poor.  However,  for  the  1 1  runs  for  which  data  for  the 
neck  were  available,  only  two  of  the  11  X  displacement  variables  were  rated  as  excellent,  with  the  others 
rated  as  good.  In  addition,  with  the  exception  of  run  LX3852,  the  magnitude  of  the  corrections  to  the 
T1  accelerometers  were  all  classified  as  excellent. 

It  may  be  noted  from  Table  D-2  that  a  less  than  excellent  rating  for  the  Marjorized  linear  and  angular 
displacement  components  was  for  the  most  part  observed  only  in  the  X  displacement  components  for  both 
the  head  and  neck,  particularly  for  the  neck.  The  tendency  of  the  X  components  for  the  head  and  the 
neck  to  be  rated  as  less  than  excellent  is  almost  certainly  the  result  of  the  previously  noted  timing 
discrepancy,  and  it  was  previously  shown  that  improved  results  can  be  obtained  by  applying  a  small 
timing  shift  to  the  sled  variables  (Figure  C-16  through  Figure  C-19). 

Whenever  a  Marjorized  variable  is  judged  to  be  poor,  an  effort  should  be  made  to  improve  the  results. 
The  possibility  of  bad  data  should  be  considered  first,  since  this  could  be  the  source  of  the  unsatisfactory 
results.  In  addition,  some  of  the  input  data  values  may  be  manipulated  to  improve  the  agreement  between 
the  Marjorized  and  the  photo-derived  variables. 

An  effort  was  made  to  improve  the  agreement  between  the  head  Marjorized  and  photo-derived 
displacement  variables  for  nm  LX3970,  since  the  head  Z-displacement  component  was  rated  as  poor. 
In  addition,  although  LX3904  did  not  have  any  components  rated  as  poor,  it  had  two  head  components 
which  were  less  than  excellent,  and  work  on  LX3904  was  therefore  warranted.  An  effort  was  also  made 
to  improve  the  Marjorized  neck  results  for  LX3959,  since  two  of  the  displacement  components  were 
classified  as  good.  For  these  runs  the  factors  which  were  applied  to  the  a  priori  variances  for  the  photo- 
derived  variables  were  manipulated  (Table  D-12).  The  poor  result  was  upgraded  to  excellent,  and  three 
of  the  five  variables  originally  categorized  as  good  were  upgraded  to  excellent.  Note,  however,  that  the 
magnitude  of  the  corrections  generally  increased,  although  not  to  an  unacceptable  level  for  the  factors 
used. 

As  previously  discussed,  a  timing  shift  was  applied  to  the  sled  acceleration  for  LX3904  and  LX3965 
in  an  effort  to  improve  the  X  displacement  components  for  the  head  and  neck.  It  may  be  seen  in  the 
bottom  plots  in  Figure  C-16  through  Figure  C-19,  which  reflect  this  timing  change,  that  a  1-msec  timing 
shift  resulted  in  much  improved  Marjorized  X  displacement  components.  Table  D-13  presents  an 
evaluation  of  the  results  for  these  timing-shift  runs.  It  may  be  seen  that  all  of  the  X  displacement 
components  were  upgraded  to  excellent. 
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Detailed  Instructions  for  (Generating  Marjorized  Variables 

Detailed  instructions  for  executing  the  Marjorization  program  on  the  HP9000  are  provided  below. 
Also  provided  are  instructions  for  gathering  the  required  data,  preparing  it  for  the  Marjorization  program, 
and  creating  the  required  data  files.  Directions  for  generating  plots  comparing  Marjorized,  photo-derived, 
and  sensor-derived  displacement  variables  are  included.  These  variables  are  defined  in  Table  D-9. 

The  flowchart  presented  in  Figure  C-1  shows  the  input  data  required  by  each  of  the  Marjorization 
modules  and  by  the  auxiliary  modules  for  run  LX3959.  It  also  indicates  the  output  from  each  module 
which  is  used  by  follow-on  modules.  Table  D-1  summarizes  the  input  and  output  data  for  the  various 
Marjorization  modules  and  lists  the  tables  that  contain  descriptions  of  the  variables  and  the  figures  that 
contain  sample  data  files  for  the  various  modules. 

Collect  Run  Data.  The  first  step  in  generating  Marjorized  variables  is  to  collect  the  data  for  the  desired 
run.  If  reduced  EZFLOW  photo  variables  for  the  run  are  available,  it  is  necessary  to  identity  the 
directory  in  which  they  are  located.  This  directory  should  be  defined  as  unit  7;  i.e.,  file  7,  located  in 
the  directory  on  the  HP  from  which  the  Marjorization  run  is  to  be  started.  File  7  should  consist  of  a 
single  entry,  the  name  of  the  directory  where  the  reduced  EZFLOW  variables  are  located.  If  the  reduced 
EZFLOW  variables  are  not  available,  it  is  necessary  to  run  EZFLOW  to  compute  them.  Sensor 
accelerometer  data  should  be  available  on  unit  8;  i.e.,  file  8  should  contain  only  the  name  of  the  directory 
where  the  accelerometer  data  is  located. 

For  example,  the  Marjorization  program  was  executed  for  run  LX3959  from  the 
/7933/users/margie/MARJORIZATION  directory  with  input  data  shown  in  Appendix  B3  through  B7  and 
Appendix  1 1 .  File  7  in  the  directory  /7933/users/margie/MARJORIZATION  contained  the  entry  (without 
quotes)  "/7963B/prod/data/esvpaper”  —  the  directory  in  which  the  reduced  EZFLOW  variables  (photo 
and  sensor)  for  LX3959  were  stored.  Unit  8  contained  the  entry  77963B/prod/data/sensor"  —  the 
directory  in  which  the  accelerometer  data  for  LX3959  were  stored. 

Reformat  Reduced  Photo  Variables.  It  is  necessary  to  reformat  the  reduced  photo  variables  to  an 
ASCII  file  to  be  read  into  the  second  Marjorization  module.  Ronnie  Wilson  wrote  the  subroutine 
getphodatf  for  this  purpose.  The  driver  program  niarjphopgm.f  reads  a  data  file,  such  as  marjphodat 
(Appendix  B8)  from  standard  input.  This  data  file  consists  of  only  two  entries:  the  unit  that  contains  the 
name  of  the  directory  in  which  the  reduced  photo  variables  are  stored  and  the  run  number.  For  example, 
for  LX3959  marjphodat  contains  the  entry  “7”  and  the  run  “LX3959,”  since  file  7  contains  the  directory 
in  which  the  reduced  photo  variables  are  stored,  i.e.,  /7963B/prod/data/esvpaper.  A  description  of  the 
variables  in  this  input  file  is  provided  in  Table  D-8. 

The  reformatted  photo  variables  are  output  to  files  with  names  constructed  with  eight  characters,  the 
first  four  of  which  are  either  “head”  or  “neck,”  and  the  last  four  characters  consist  of  the  four-digit  run 
number.  For  example,  for  run  LX3959  the  reformatted  data  files  are  head3959  (Appendix  B9)  and 
neck3959. 

The  executable  file,  xqtmarjpho,  which  reformats  the  photo  variables,  is  located  in  the 
/7933/users/ronnie/MARJORIZATION  directory.  The  program  can  be  executed  with  the  xqtpho  script 
shown  in  Appendix  B2,  or  by  entering  the  following  command,  assuming  that  the  data  file  marjphodat 
contains  the  appropriate  input  data: 


7 


_ Naval  Biodynamics  Laboratory  Software  Documentation _ 

/7933/users/ronnie/MARJORIZATION/xqtmarjpho<marjphodat 

Marjorization  Program.  The  Marjorization  program  can  be  executed  for  any  sled  run 
using  a  customized  version  of  the  shell  script  xqtmarj  (Appendix  Bl).  Before  xqtmarj  is  wecuted 
owever,  it  is  necessary  to  prepare  the  individual  scripts  and  data  files  summarized  in  Table  D-1  and 

Sample  scripts  and  user-defined  input  data  files  used  to  execute  the  Marjorization  program  for  LX3959 
HPO^n”  I"  ‘""“S'’  “‘I  Bll.  These  files  are  currZy”L  ™ 

fndTnni"  T  1 — /7933/users/margie/MARJORIZATION.  If  the  data  specified  above  are  available 
a  H  f  baseline  input  data  are  to  be  used  for  the  run,  then  the  only  changes  required  to  the  scripts 

mes«'  Zaj  tr"f  four-digit\led  run  number'in 

2an  nm  and  So/  tt  "’’^^omical  mounts  for  the  head  and  T1  are  other 

than  1101  and  2201  it  is  also  necessary  to  provide  perturbation  files  for  the  anatomical  mounts  used 

similar  to  the  files  shown  in  Appendix  B4.  (See  Table  D-3  for  a  definition  of  the  Za  varTawS  S 

ming  convention  for  these  files  is  “d”  concatenated  with  the  four-digit  anatomical  mount  number 

^  scripts  and  data  files  have  been  prepared,  the  Marjorization  program  may  be  executed  bv 

entering  the  following  at  the  “$”  prompt  on  the  HP9000:  ^ 

xqtmarj 

The  script  xgtmarj  provides  for  Ihe  execution  of  all  three  primary  Marjorization  modules  Iterative 

r„e Tf  °ttained  wi,h  r  cLm^d 

for  i  in  corrl  corr2 

If  three  iterations  are  desired,  then  the  user  should  append  a  space  and  “corr3”  to  the  above  line.  Anv 
desired  number  of  iterations  can  be  obtained  similarly.  ^ 

The  sequence  of  operations  specified  in  the  script  xqtmarj  is  as  follows: 

.  nor  to  the  start  of  the  iteration  loop,  the  script  copies  the  file  corrzero  (containin<T  all  zero 
wrrections)  to  the  file  corr3959  to  initialize  the  correction  file  for  the  first  iteration 
.  script  then  initiates  a  loop  in  which  script  rwil  (Appendix  B2)  executes  Marjorization 
Module  1  using  the  corrections  to  the  sensor  parameters  found  in  file  corriP5P  to  compute 
reference  values  for  the  sensor  parameters  to  which  perturbations  are  applied.  The  output  from 
Marjorization  Module  1  is  the  file  out3959,  which  contains  perturbed  and  unperturbed  sensor 
to^'^irid  (Appendix  Bl)  then  copies  the  correction  file 

),  which  executes  Marjorization  Module  2.  Corrections  from  oldc3959  are  applied  to  the 

the  reference  values  about  which  the  constrained  least-squares 
solution  IS  linearized  to  compute  corrections  to  the  sensor  parameters  for  that  iteration 

f  for IrnexTSTon.  Module 

(XpetduT)' thtr"'’’  T"  M  ">0  "-“'or  scrip,  xgmar/  runs  scrip,  runS 

(Appendix  B2),  which  executes  Marjorization  Module  3.  Module  3  reads  the  corrections 

contained  in  corr3959,  applies  them  to  the  sensor  parameters,  and  generates  corrected  sensor 
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variables,  which  are  defined  as  Marjorized  variables. 

Detailed  information  on  the  execution  of  the  individual  Marjorization  modules  is  presented  below. 

Execution  of  Marjorization  Module  1.  Execution  of  Marjorization  Module  I  may  be  accomplished 
by  executing  the  script  runl  shown  in  Appendix  B2.  This  runs  the  executable  program  xqtmarjptli 
located  in  the  directory  /7933/users/ronnie/MARJORIZATION.  This  program  requires  File  7  to  contain 
the  name  of  the  directory  in  which  the  EZFLOW  reduced  photo  variables  (and  sensor  variables)  for  the 
given  run  are  stored,  e.g.,  /7963B/prod/data/esvpaper  for  LX3959  data.  File  8  must  contain  the  name 
of  the  directory  in  which  the  accelerometer  data  are  stored,  e.g.,  /7963B/prod/data/sensor  for  LX3959 
sensor  data. 

Files  containing  deltas  to  be  applied  to  the  head  and  neck  mount  accelerometer  direction  cosines, 
sensitivities,  and  initial  Euler  angles  to  generate  perturbed  values  must  be  available  in  the  current 
directory.  A  file  must  be  provided  for  each  mount,  and  the  file  name  must  be  constructed  as  “d” 
concatenated  with  the  four-digit  anatomical  mount  number.  For  example,  if  the  head  and  neck  anatomical 
mounts  are  1101  and  2201  respectively,  then  the  deltas  must  be  stored  in  files  named  dllOl  and  d2201 
(Appendix  B4).  Data  required  from  standard  input  are  described  in  Table  D-4,  and  sample  data  are 
shown  in  data  file  datal  in  Appendix  B5. 

In  order  to  permit  iterative  solutions  with  linearization  around  the  solution  from  the  previous  iteration, 
the  program  reads  in  corrections  to  be  applied  to  the  accelerometer  parameters  for  each  iteration.  On  the 
first  iteration,  the  program  must  read  all  zero  corrections.  Hence,  a  data  file  in  which  all  corrections  are 
set  to  0  must  be  provided,  e.g.,  corrzero  (Appendix  B1 1).  On  line  2  of  xqtmarj,  this  data  file  with  zero 
corrections  is  copied  to  the  correction  file  to  be  used  for  that  iteration,  e.g.,  corr3959.  The  name  of  this 
correction  file  is  the  concatenation  of  “corr”  and  the  four-digit  run  number,  e.g.,  corr3959  for  run 
LX3959.  On  successive  iterations  the  correction  file  contains  computed  corrections  output  by 
Marjorization  Module  2  on  the  previous  iteration. 

Output  from  Marjorization  module  1  is  a  file  containing  the  perturbed  and  unperturbed  values  for  the 
sensor-derived  linear  displacements  and  quaternions.  This  data  file  contains  unperturbed  and  perturbed 
sensor-derived  variables  obtained  by  applying  the  perturbation  values  shown  in  Appendix  B4  and 
described  in  Table  D-3.  -This  file  contains  identifying  information  such  as  run  number,  data  type  (head 
or  neck),  number  of  data  points,  time  increment,  and  the  magnitude  of  the  perturbation  applied.  For  each 
displacement  variable  it  contains  both  the  unperturbed  value,  and  the  perturbed  value  for  each 
perturbation  applied  to  a  sensor  parameter.  The  file  name  is  the  concatenation  of  “out”  and  the  four-digit 
run  number,  e.g.,  out3959.  The  program  assumes  that  this  file  does  not  exist.  If  it  does,  it  is  necessary 
to  delete  it. 

Sample  script  rml  in  Appendix  B2  may  be  used  to  execute  this  module  for  LX3959.  Note  that  the 
master  script  xqtmarj  (Appendix  Bl)  contains  a  command  to  delete  the  fi\&  out3959  if  it  exists.  Script 
rml  contains  the  command  to  execute  Marjorization  module  1  with  data  for  LX3959: 

/7933/users/ronnie/MARJORIZATION/xqtmarjptl<datal 

In  addition  to  the  writing  data  to  out3959,  reduced  EZFLOW  unperturbed  sensor  variables  named  with 
lower-case  letters  are  written  to  unit  7  by  this  module. 

Execution  of  Marjorization  Module  2.  Execution  of  Marjorization  Module  2  may  be  accomplished 
by  executing  the  script  rm2  (Appendix  B2).  Marjorization  Module  2  requires  an  input  file  whose  name 
is  the  concatenation  of  “oldc”  and  the  four-digit  run  number,  e.g.,  oldc3959,  which  contains  corrections 
from  the  previous  iteration.  (For  the  first  iteration,  of  course,  the  corrections  would  be  all  zeroes).  The 
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format  for  this  file  is  shown  in  Appendix  BIO.  Module  2  also  requires  the  Module  1  output  file  “out” 
plus  the  four-digit  run  identification  number  (“runid”),  as  well  as  reduced  head  and  neck  photo  variables, 
which  must  be  stored  in  files  named  “head”  plus  the  four-digit  runid  and  “neck”  plus  the  four-digit  runid. 
For  run  LX3959  these  files  are  out3959,  hecid3959,  and  neck3959.  The  program  also  requires  that  there 
not  be  an  existing  correction  file,  such  as  corr3959,  since  this  file  is  produced  by  the  program. 

Master  script  xqtmarj  (Appendix  Bl)  contains  commands  to  copy  the  existing  correction  file  from  the 
previous  iteration,  corr3959,  to  the  input  file  oldc3959,  and  to  delete  the  file  corr3959  as  required  by  the 
program.  Script  run2  contains  the  command  to  run  the  executable  program  xqtmarjpt2i,  which  is  located 
in  the  /7933/users/ronnie/MARJORIZATION  directory.  Correction  values  read  from  oldc3959  are 
applied  to  the  sensor  parameters,  which  become  the  reference  values  about  which  the  constrained  least- 
squares  solution  is  linearized  to  compute  corrections  to  the  sensor  parameters  for  that  iteration.  The 
program  reads  the  f\\tfllein  (Appendix  B6)  from  standard  input.  A  description  of  the  input  variables  in 
fllein  is  provided  in  Table  D-5.  Filein  contains  the  name  of  the  data  file  that  contains  additional  run  data, 
such  as  data2  (Appendix  B6),  as  well  as  factors  to  be  applied  to  the  denominators  of  the  variances  for 
the  a  priori  photo-derived  displacement  variables.  A  description  of  the  input  data  variables  in  data2  is 
provided  in  Table  D-6.  Note  that  this  method  of  entering  the  data  was  used  to  facilitate  use  of  the 
debugger  without  having  to  enter  the  Marjorization  data  from  the  keyboard. 

Output  from  Marjorization  Module  2  is  a  file  named  “corr”  plus  the  four-digit  runid,  e.g.,  corr3959, 
which  contains  the  cumulative  corrections  to  the  direction  cosines,  sensitivities,  and  initial  Euler  angles 
for  both  the  head  and  neck,  which  are  computed  by  Marjorization  Module  2.  These  corrections  are 
required  by  Marjorization  Module  3  and  are  also  read  by  Marjorization  Module  1  when  iterative  solutions 
are  obtained. 

Marjorization  Module  2  may  be  executed  without  using  script  run2  by  entering  the  following  command 
after  the  appropriate  data  files  have  been  provided  as  described  above: 

/7933/users/ronnie/MARJORIZATION/xqtmar jpt2i<filein 

Execution  of  Marjorization  Module  3.  Execution  of  Marjorization  Module  3  may  be  accomplished 
by  executing  the  script  run3  (Appendix  B2).  This  script  runs  the  executable  program  xqtmarjpt3,  which 
is  located  in  directory  /7933/users/ronnie/MARJORIZATlON.  This  program  requires  three  lines  of  input 
from  the  standard  input  device,  as  shown  in  file  dataS  in  Appendix  B7.  A  description  of  data  variables 
is  provided  in  Table  D-7.  The  program  also  requires  a  file  with  corrections,  such  as  corr3959  (Appendix 
BIO),  which  is  output  by  Marjorization  Module  2,  to  be  applied  to  the  accelerometer  parameters.  It  also 
requires  access  to  reduced  photo  variables  on  unit  7  and  accelerometer  data  on  unit  8,  as  discussed  for 
Marjorization  Module  1. 

Output  from  this  module  consists  of  Marjorized  variables,  which  are  corrected  sensor-derived  variables 
named  with  upper-case  letters,  and  with  a  final  letter  “C.”  The  Marjorized  variables  are  stored  in  the 
directory  specified  in  file  7. 

The  following  command  may  be  used  to  execute  this  module,  if  script  rm3  is  not  used,  once  the 
required  data  files  have  been  provided: 

/7933/users/ronnie/MARJORIZATION/xgtmar jpt3<data3 
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Plotted  Output.  Plots  of  Marjorized  variables  may  be  obtained  using  the  usual  Data  Systems  plot 
programs.  Sample  plotting  elements  as  well  as  shell  scripts  to  automatically  generate  plots  are  stored  in 
the  directory  /7933/users/margie/MARJORIZATION/PLOTS. 

Standard  plotted  output  for  each  Marjorization  run  is  a  set  of  12  plots  comparing  Marjorized  linear  and 
angular  displacement  variables  (file  name  ending  with  the  letter  “C”),  photo-derived  variables  (final  letter 
“P”),  and  sensor-derived  variables  (final  letter  “S”).  Definitions  of  the  plotted  variables  are  provided 
in  Table  D-9.  Figure  C-2  through  Figure  C-13  are  a  standard  set  of  plots  for  LX3959. 

These  plots  may  be  generated  using  a  customized  version  of  the  plot  specification  file  pltmarj,  which 
is  stored  on  the  HP9000  in  the  directory  /7933/users/margie/MARJORIZATION/PLOTS.  Customization 
of  this  file  consists  of  replacing  the  four-digit  run  number  in  the  existing  file  with  the  desired  run  number 
and  changing  the  title  on  line  39  as  appropriate.  In  addition,  the  units  on  which  the  variables  to  plot  are 
located  must  be  provided  on  lines  9  and  10  of  the  file  pltmarj.  For  example,  for  LX3959  line  9  contained 
the  entry  “1”  and  file  1  contained  the  entry  "/7963B/prod/data/esvpaper,''  the  directory  in  which  the 
Marjorized,  photo-derived,  and  sensor-derived  variables  were  stored. 

Hie  script  autoplt  shown  in  Appendix  B2  may  be  used  to  automatically  generate  these  plots  to  the 
screen  of  an  HP2627  terminal;  it  is  set  up  to  read  the  post-processor  file  pop2627  in  the  above  directory. 
To  generate  hardcopy  plots,  the  autoplt  script  may  be  modified  to  read  the  post-processor  file pop7550; 
then  at  the  “$”  prompt  this  command  may  be  entered: 

$plot7550 

Alternately,  after  generating  screen  plots,  hardcopy  plots  may  be  executed  by  entering 

$xgt_pop<pop7 5 5 0 


and  then  entering 


$plot7550 


Summary  and  Conclusions 

Marjorized  variables  suitable  for  mathematical  modeling  may  be  obtained  by  using  the  procedures 
described  in  this  document.  Detailed  instructions  for  generating  Marjorized  variables  have  been  provided 
above,  and  nominal  values  for  all  input  data  appear  in  the  accompanying  tables  and  appendices.  The 
baseline  data  values  described  in  this  paper  are  appropriate  for  production-type  runs.  Marjorized  results 
may  be  evaluated  using  the  criteria  specified  in  Table  D-2.  If  the  evaluation  of  Marjorized  variables 
generated  with  the  baseline  data  values  results  in  less  than  the  desired  agreement  between  the  photo- 
derived  and  Marjorized  displacement  variables,  it  would  be  appropriate  to  first  investigate  the  possibility 
that  there  is  a  problem  with  the  data,  and  then  to  manipulate  input  data  as  previously  described  to 
improve  this  agreement. 

Several  extensions  to  the  Marjorization  program  are  desirable  for  future  implementation.  The 
capability  of  computing  corrections  to  the  angular  velocity  sensors  would  be  useful  in  Marjorizing  the 
more  recent  sled  runs,  which  employ  angular  velocity  sensors.  In  addition,  the  capability  to  input 
variable  time  shifts  for  the  sled  acceleration  trace  should  be  implemented  as  an  option. 
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Implementation  of  the  Maijorization  Program 


Assumptions: 


Marjorized  Variables: 


Corrections: 


Constraints: 


Appendix  A 

Summary  of  the  Marjorization  Process 


(1)  EZFLOW  photo-derived  displacements  are  the  best  source  of  displacement 
variables. 

(2)  EZFLOW  sensor-derived  accelerations  are  the  best  source  of  acceleration 
variables. 

(3)  Small  corrections  to  various  accelerometer  constants  can  be  determined;  these 
corrections: 

(a)  Produce  very  little  change  in  the  acceleration  variables 

(b)  Can  produce  significant  changes  in  the  computed  displacements 

(c)  Bring  the  corrected  sensor-derived  displacements  into  agreement  with  the 
photo-derived  displacements 

The  corrected  accelerometer-derived  kinematic  variables  obtained  by  applying  the 
computed  corrections  to  the  appropriate  parameters  in  the  EZFLOW  sensor 
program.  These  provide  a  consistent  set  of  kinematic  variables  from  the 
acceleration  to  the  displacement  level. 

The  following  corrections  are  intended  to  bring  the  accelerometer-derived  linear 
displacements  and  angular  orientations  into  agreement  with  the  corresponding 
photo-derived  variables: 

(1)  Accelerometer  direction  cosines 

(2)  Accelerometer  sensitivities 

(3)  Initial  Euler  angles 

The  algorithm  includes  constraint  terms  intended  to  deter  excessively  large 
computed  correction  components. 
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Appendix  B 

Sample  Scripts  and  Data  Files 


Bl.  Shell  Script  xqtmaij  to  Execute  Marjorization  Program  for  LX3959  with  Two  Iterations 

(See  Appendix  B2  for  the  individual  module  scripts.) 


#xqtmar j 

cp  corrzero  corr3959  #initialize  correction  components  to  0. 
for  i  in  corrl  corr2 
do 

echo  correction  no.  $i 

rm  out3959  /delete  old  file  -  this  is  output  by  Module  1 

runl>/7914/atbmodel/outl 

cp  corr3959  oldc3959  /copy  old  correction  file  to  oldc3959 

rm  corr3959  /delete  old  file  ~  this  is  output  by  Module  2 

run2>/7914/atbmodel/out2 

cp  corr3959  $i  /store  corrections  in  files  corrl,  etc. 

done 

run3>out3 

/print  data  files  and  final  output 

Ip  xqtmarj  runl  run2  datal  data2  filein  corrl  run3  data3  corr3959  out3 


B2.  UNIX  Shell  Scripts  to  Execute  Individual  Marjorization  Modules 

script  runl  —  Marjorization  Module  1: 

#execute  Module  1  with  data  file  datal 
/7933/users/ronnie/MARJORIZATION/xqtmarjptli<datal 


script  nin2  —  Marjorization  Module  2: 

#execute  Module  2  with  data  file  filein 
/7933/users/ronnie/MARJORIZATION/xqtmarjpt2i<filein 


script  runs  —  Marjorization  module  3: 

#execute  Module  3  with  data  file  data3 
/7933/users/ronnie/MARJORIZATION/xqtmarjpt3<data3 
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script  xqtpho  —  reformat  photo-derived  variables  required  by  module  2: 

#execute  photo-reformatting  module  with  data  file  marjphodat 
/ 79  3  3 /user s /ronnie /MARJORI ZATION/xqtmar j  pho<mar j  phodat 


script  autoplt  —  generate  comparison  plots  of  photo-derived,  accelerometer-derived,  and  Marjorized 
displacement  variables: 

$xqt_specs  <  pltmarj  >  out  #plot  specs  from  file  pltmarj 

tail  out 

$xqft_plot 

mv  popfil.dat  meta.dat 

$xgt_pop  <  pop2627  #display  plots  on  screen  of  HP2627A  terminal 


B3.  Baseline  Input  Data  Files  for  Marjorization  Module  2 

(See  Table  D-5  for  definitions  of  data  entries  in  filein  and  data2,) 

Data  file  filein: 

data2 

1.0  1.0  1.0 

3.0  3.0  3.0 


Data  file  data2: 

LX3959  0  1  1.0000  1.0000 
0.0700  0.1000  0.1800  0.2300  0.2500  0.2900  0.3090 

0.0700  0.1000  0.1800  0.2300  0.2500  0.2900  0.3090 

1  2  2  2  2  2  1 
1  2  2  2  2  2  1 


10.0  10.0  10.0 

3.0  3.0  3.0 
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B4.  Sample  Data  Files  for  Input  to  Marjorization  Module  1  —  Perturbation  Values  for 
Accelerometer  Parameters 

(See  Table  D-3  for  definitions  of  data  in  dllOl  and  d2201.) 

Data  file  dllOl: 

amnt  1101  DELTAS  FOR  MARJORIZATION 
0.0087178 
0.0087178 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0087178 
0.0087178 
0.0060000 
0.0087178 
0.0087178 
0.0060000 
0.0061644 
0.0061644 
..  0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0435890 
0.0435890 
0.0435890 
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Data  file  d2201: 

amnt  2201  DELTAS  FOR  MARJORIZATION 
0.0087178 
0.0087178 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0087178 
0.0087178 
0.0060000 
0.0087178 
0.0087178 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
0.0061644 
0.0061644 
0.0060000 
..  0.0061644 

0.0061644 
0.0060000 
0.0435890 
0.0435890 
0.0435890 


B5.  Sample  Input  Data  File  for  Marjorization  Module  1 

(See  Table  D-4  for  definitions  of  the  data  entries.) 

Data  file  datal: 

1  0 

LX3959 

EOFEOF 
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B6.  Original  Marjorization  Module  2  Input  Data  Files  for  LX3959  —  Used  for  Plots  in 
Figure  C-2  through  Figure  C-13. 

(See  Table  D-5  for  definitions  of  data  entries  in  filein  and  data2.) 

Data  file  filein'. 
data2 

1.0  1.0  1.0  10.0  10.0  10.0 

5.0  5.0  5.0  5.0  5.0  5.0 


Data  file  data!-. 

LX3959  0  1  3.0000  3.0000 
0.1200  0.1600  0.2000  0.2600  0.2850  0.2900  0.2990 

0.1200  0.1600  0.2000  0.2600  0.2850  0.2900  0.2990 

12  3  2  1  2  2  2 

12  3  2  1  2  2  2 


B7.  Sample  Input  Data  File  for  Marjorization  Module  3 
(See  Table  D-6  for  definitions  of  data  entries  in  dataS.) 

Data  file  dataS: 

1  0 

LX3959 

EOFEOF 


B8.  Sample  Input  Data  for  Photo-Reformatting  Marjorization  Module 

(See  Table  D-7  for  definitions  of  data  entries  in  marjphodat.) 

Data  file  marjphodat. 

7  LX3959 
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B9.  Sample  Output  from  Photo-Reformatting  Module 
(See  Table  D-10  for  a  description  of  the  output  data  variables.) 

Initial  records  from  output  data  file  head3956: 


time 

.OOOOOOOOE  +  00 

.20999999E-02 

.4i99ggg9E-02 

.62000002E-02 

.82999999E-02 

.10400000E-01 

.12600000E-01 

.14600000E-01 

.16600000E-01 

.18700000E-01 

.20800000E-01 

.22900000 E-01 

.24900001  E-01 

.27000001  E-01 

.29 10000  IE-01 

,31 199999E-01 

.33199999E-01 

.36300002E-01 

.37400000E-01 

.39600002E-01 

.41600000E-01 

.43600000E-01 

.46699999E-01 

.47800001  E-01 

,49899999E-01 

.61899999E-01 

.6400000  IE-01 

.661 00000 E-01 

.68200002E-01 

.60199998E-01 

.e2300000E-01 

.64400002E-01 

.66600001  E-01 

.68499997E-01 

.70600003E-01 

.72700001  E-01 

.74800000E-01 

.76800004E-01 

.78900002E-01 


daxsop 

-.12818400E  +  01 
-.12816900E  +  01 
-.12816600E-I-01 
-.12814000E  +  01 
-.12812300E  +  01 
-.12810900E-»-01 
-.12810100E-I-01 
-.12809900E-h01 
-.12810000E-H01 
-.12809900E  +  01 
-.12809800E  +  01 
-.12809600E-1-01 
-.12808900E-1-01 
-.12808200E  +  01 
-.12807300E  +  01 
-.12806300E-H01 
-.12805400E  +  01 
-.12806001E  +  01 
-.12804700E-J-01 
-.1 2803800  E-h  01 
-.12802401E  +  01 
-.12801000E  +  01 
-.12799600E  +  01 
-.12797700E-h01 
-.12794B01E  +  01 
-.12790200E  +  01 
-.12784801E  +  01 
-.12778000E  +  01 
-.12769099E  +  01 
-.12767600E  +  01 
-.12742200E  +  01 
-.12722300E  +  01 
-.12697600E  +  01 
-.12667700E  +  01 
-.12631600E  +  01 
-.12B89200E  +  01 
-.12B41800E-^01 
-.12490100E  +  01 
-.12434700E-(-01 


daysop 

-.61871698E-02 
-.62933699E-02 
-.64227101 E-02 
-.66447602E-02 
-.66269902E-02 
-.66442299E-02 
-.662BB000E-02 
-.6676690  IE-02 
-.6604630  IE-02 
-.64474801  E-02 
-.64249299E-02 
-.64346398E-02 
-.64749299E-02 
-.6B434698E-02 
-.66866 100E-02 
-.6569730  IE-02 
-.66266601  E-02 
-.64823702E-02 
-.6432820  IE-02 
-.63623800E-02 
-.62804399E-02 
.-.62606699E-02 
-.63040499E-02 
-.642 1680  IE-02 
-.66201 302E-02 
-.66699699E-02 
-.66446601  E-02 
-.67893001  E-02 
-.69640602E-02 
-.71271299E-02 
-.7264360  IE-02 
-.72960602E-02 
-.72173099E-02 
-.71332902E-02 
-.71297898E-02 
-.73638898E-02 
-.78679603E-02 
-.8B377796E-02 
-.92312200E-02 


dazsop 

.16472900E  +  01 
.16474700E  +  01 
.1B476000E-h01 
.16477300E  +  01 
.16478600E  +  01 
.16480100E-t-01 
.16481400E  +  01 
.16481600E  +  01 
.164814006  +  01 
.16481600E  +  01 
.16482600E  +  01 
.164840006  +  01 
.1B48B300E  +  01 
.16486400E  +  01 
.16487400E  +  01 
.16488000E  +  01 
.16488300E  +  01 
.16488300E  +  01 
.16488300E  +  01 
.16488300E  +  01 
.1B488601E  +  01 
.16489200E  +  01 
.16490100E  +  01 
.16491000E  +  01 
.16491700E  +  01 
.16492100E  +  01 
,164922006  +  01 
.16492400E  +  01 
.1B493000E  +  01 
.16494000E  +  01 
.16494800E  +  01 
.16496200E  +  01 
.16494800E  +  01 
.16492899E  +  01 
.16489399E  +  01 
.16484101E  +  01 
.16476800E  +  01 
.16467300E  +  01 
.16466000E  +  01 


4hoo1p 

.31407000E-02 
.31 764901  E-02 
.32420601  E-02 
.336661 99E-02 
.34829699E-02 
.36691 799E-02 
.36641400E-02 
.33876901  E-02 
.32117399E-02 
,30743701  E-02 
.30203301  E-02 
.30682 100E-02 
.31604299E-02 
.32834 10  IE-02 
.33284600E-02 
.33163701E-02 
.32704901  E-02 
.31415699E-02 
.30280000E-02 
.29947199E-02 
.29461  300E-02 
.296103996-02 
.30764299E-02 
.33661 101E-02 
.37171999E-02 
.40328498E-02 
.42569200E-02 
.44289 199E-02 
.47370102E-02 
.B2094699E-02 
.669143986-02 
.60308501  E-02 
.61889901E-02 
.61347801E-02 
.688B7398E-02 
,667461996-02 
.B3333300E-02 
.60181202E-02 
.44447999E-02 


4hoo2p 

.460327996-01 
.46814399E-01 
.46612399E-01 
.47111899E-01 
.476106006-01 
.48112199E-01 
.48566101  E-01 
.48746200E-01 
.48864899E-01 
.49071699E-01 
.494666006-01 
,499443006-01 
.60376698E-01 
.60768901  E-01 
.B1139600E-01 
.61481 600E-01 
.61761200E-01 
.61876200E-01 
.619820006-01 
.621233016-01 
.62233100E-01 
.62438200E-01 
.62690301  E-01 
.62961 101  E-01 
.63196398E-01 
.63384699E-01 
.636690986-01 
.63767700E-01 
.64104000E-01 
,646374996-01 
.66267200E-01 
.66903602E-01 
.66409702E-01 
.66601  800E-01 
.66669900E-01 
.664663026-01 
.663413996-01 
.66299899E-01 
.66401901E-01 


4hoo3p 

-.60249101  E-01 
-.60080 100E-01 
-.698989996-01 
-.69710801  E-01 
-.69601 299E-01 
-.69664600E-01 
-.69780400E-01 
-.69967801  E-01 
-.60246302E-01 
-.60646899E-01 
-.60688499E-01 
-.60641401  E-01 
-.60644 100E-01 
-.60396600E-01 
-.603483996-01 
-.60614499E-01 
-.606492016-01 
-.60647000E-01 
-.60719401E-01 
-.60986200E-01 
-.61304699E-01 
•.61341699E-01 
-.61173800E-01 
-.60996499E-01 
-.61008200E-01 
-.61169998E-01 
-.61146000E-01 
-.60949899E-01 
-..60766999E-01 
-.60704999E-01 
-.60934 100E-01 
-.61694799E-01 
-.62691 197E-01 
-.640382026-01 
-.66383904E-01 
-.66281497E-01 
-.66467698E-01 
-.6601 1101  E-01 
-.6B122098E-01 


4hoo4p 

.99716799E  +  00 
.99713397E  +  00 
.9971 1 400  E  + 00 
.99709803E  +  00 
,997080986  +  00 
.99706303E  +  00 
.997024006  +  00 
.99700397E  +  00 
.996981986  +  00 
.99696402E  +  00 
.99692601E  +  00 
.996904026  +  00 
.99688798E  +  00 
.99687701 E  + 00 
.99686098E  +  00 
.99683303E  +  00 
.99681097E  +  00 
.99680601  E  +  00 
.99679601E  +  00 
.99677199E  +  00 
.99674702E  +  00 
.99673402E  +  00 
.99673003E  +  00 
.99672698E  +  00 
.99671 102E  + 00 
.99668998E  +  00 
.996680976  +  00 
.99668097E  +  00 
.99667299E  +  00 
.99664497E  +  00 
.99669300E  +  00 
.99661402E  +  00 
.99641 697E +  00 
.99631 900E  + 00 
.99623601  E  +  00 
.99618600E  +  00 
.99618101E  +  00 
.99621 600E  + 00 
.99627000E  +  00 
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BIO.  Computed  Corrections  for  LX3959,  Iteration  2 

Data  file  corr3959  must  be  initialized  to  all  zero  components  prior  to  executing  Marjorization  Module 
1.  This  is  accomplished  by  copying  the  file  corrzero  (Appendix  Bll)  into  the  file  corr3959.  The  script 
xqtmarj  in  Appendix  B1  accomplishes  this  by  the  command  on  line  2  of  this  script.  In  this  way 
Marjorization  Module  1  does  not  apply  any  corrections  on  the  first  iteration.  The  correction  file  corr3959 
must  be  copied  to  the  file  oldc3959  (representing  corrections  from  the  previous  iteration),  and  corr3959 
must  be  deleted,  since  this  file  is  output  by  Marjorization  Module  2.  This  is  accomplished  automatically 
by  the  commands  on  lines  8  and  9  of  sciv^t  xqtmarj  .  Marjorization  Module  2  then  computes  corrections 
for  the  current  iteration  which  are  stored  in  corr3959.  This  file  is  input  to  Marjorization  Module  1  for 
successive  iterations.  The  file  corr3959  output  by  Marjorization  Module  2  on  the  final  iteration  is  then 
read  and  applied  by  Marjorization  Module  3.  (See  page  9  for  definitions  of  the  data  variables.) 


Data  file  corr3959  (also  oldc3959  for  the  next  iteration): 


LX3969  head  30 
.2429431 8  E-02 
.38662687E-03 
-.3656901 9E-01 
-,76414074E-02 
LX3969  neck  21 
.61606148E-02 
-.24462903E-02 
.46464081  E-02 


-.31496771  E-01 
-.21047810E-01 
-.32646082E-01 
-.11969826E-01 

-.14609219E-01 

-.90931068E-02 

-.62716762E-02 


-.327981  97E-01 
-.1 64961 04E-01 
.64189600E-02 
.16626632E-02 

.167367E7E-01 

-.27408332E-01 

.1184B039E-01 


-.13303082E-01 
.80997 190E-02 
.66972673E-02 
.1 1073916E-01 

-.31609792E-01 

.24642649E-01 

-.12629074E-01 


-.24283 137E-02 
-.11127490E-01 
.63601  730E-02 
-.17082024E-01 

.26931331E-01 

.43903468E-02 

-.64224276E-02 


.13662002E-01 
.67674682E-01 
.70709427E-03 
-.841 16291  E-02 

-.46282622E-02 

.48344713E-02 


.69324960E-02 

-.36146398E-01 

.13681426E-01 

-.6684261  6E-02 
-.34637378E-02 


.10798099E-02 

-.13671342E-01 

.40426627E-02 

-.29034908  E-02 
.10331963E-01 


Bll.  All-Zero  Corrections  for  LX3959  —  Initial  Corrections  Read  by  Marjorization 
Module  1 

(See  page  9  for  definitions  of  the  data  variables.) 


Data  file  corrzero. 


LX3959  head  30 
.OOOOOOOOE-fOO 
.OOOOOOOOE  +  00 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  00 
LX3969  neck  21 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 


.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 

.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 


.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 

.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 


.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 

.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 


.OOOOOOOOOOE  +  OO 
.OOOOOOOOOOE+00 
.OOOOOOOOOOE+OO 
.OOOOOOOOOOE  +  OO 

.OOOOOOOOOOE  +  OO 
.OOOOOOOOOOE  +  OO 
.OOOOOOOOOOE+OO 


.OOOOOOOOOOE  +  OO 
.OOOOOOOOOOE  +  OO 
.OOOOOOOOOOE  +  OO 
.OOOOOOOOOOE  +  OO 

.OOOOOOOOOOE  +  OO 
.OOOOOOOOOOE  +  OO 


.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 

.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 


.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 

.OOOOOOOOE  +  OO 
.OOOOOOOOE  +  OO 
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Appendix  C 
Figures 


Figure  C-1 

Marjorization  Program  Flowchart  for  LX3959 


Naval  Biodynamics  Laboratory  Software  Documentation 


Figure  C-2 

Comparison  of  Photo-Derived  (P),  Marjorized  (C)  and 
Sensor-Derived  (S)  Head  Anatomical  X  Displacements 


Figure  C-3 

Comparison  of  Photo-Derived  (P),  Marjorized  (C),  and 
Sensor-Derived  (S)  Head  Anatomical  Y  Displacements 
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Comparison  of  Photo-Derived  (P),  Marjorized  (C),  and 
Sensof'Derived  (S)  Head  Anatomical  Z  Displacements 


Figure  C-5 

Comparison  of  Photo-Derived  (P),  Marjorized  (C), 
and  Sensor-Derived  (S)  Rotations  About  Head 
Anatomical  X  Axis 
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F-3. OBITER  2,PS16MR  D/1,P/10 


Legend 

LX3959  PHB02P  HOP  1T!7  ]^. 7£n 

LX3959,.PmQ2C..,Hp012.7„J.1,76g.. 
LX_3959.  PHB02S  .  ,Hgg.l  27_  H .  76^ 


Figure  C-6 

amparison  of  Photo-Derived  (P),  Marjorized  (C),  and 
Sensor-Derived  (S)  Rotations  About  Carried  Head 
Anatomical  Y  Axis 


r-3.0,lTER  2,PSIGMfl  0/1, P/10 


Legend 

Lxsgsg  PHcnsp  Hnni?? 

14.76a 

LX3959  PHC03C 

HOC  127 

LX3959  PHC03S 

H00127 

14.76a 

_ — -  : - 1 

Figure  C-7 

Comparison  of  Photo-Derived  (P),  Marjorized  (C),  and 
Sensor-Derived  (S)  Rotations  About  Carried  Head 
Anatomical  Z  Axis 
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Figure  C-8 

Comparison  of  Photo-Derived  (P),  Marjorized  (C), 
and  Sensor-Derived  (S)  T1  Anatomical  X  Displacements 


F“3.0,]TER  2,PSIGMn/5 
Legend 


.LXm9  .PMSgC„.  HOO 1 27  . 76g .. 
LX3959  DNYSgS H9Q.1,27  _  A1^769. 


% 

— 

TlMt  ( SECONDS  1 


Figure  C-9 

Comparison  of  Photo-Derived  (P),  Marjorized  (C), 
and  Sensor-Derived  (S)  T1  Anatomical  Y  Displacements 
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F-3.0,ITER  2,PSIGMR/5 


Legend 

_LX3959  DNZSOP  HOP  127  14.7K«] 

.i<X3m.QNZSgC., .  . .  J.4, 76g .. 

-LX3959  D.NZSQS .  _H00,1 27,,,.  .H,,76g. 


Figure  C-10 

Comparison  of  Photo-Derived  (P),  Marjorized  (C), 
and  Sensor-Derived  (S)  T1  Anatomical  Z  Displacements 


r-3.0,IIER  2,P5IGnR/5 


Legend 

LX3959  PNRnXP 

HD0127 

l‘!.76a 

LX3959  PNROXC 

HOOl 27 

.LX.3959  P.NROXS. 

H00I27 

,  ,ll--.76g,. 

Figure  C-1 1 

Comparison  of  Photo-Derived  (P),  Marjorized  (C), 
and  Sensor-Derived  (S)  Rotations  About  T1 
Anatomical  X  Axis 
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Figure  C-1  2 

Comparison  of  Photo-Derived  (P),  Marjorized  (C),  and 
Sensor-Derived  (S)  Rotations  About  Carried  T1 
Anatomical  Y  Axis 


Figure  C-1  3 

Comparison  of  Photo-Derived  (P),  Marjorized  (C),  and 
Sensor-Derived  (S)  Rotations  About  Carried  T1 
Anatomical  Z  Axis 
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-3.0, ITER  2,P51GMR  D/1, 

P/100 

Legend 

LX3965  ORXSOP 

H00135 

M.e-Ig 

LX39B5  DRXSOC 

H00135 

...14,649.. 

..1X3965.  0R.X5Q5... 

....H00.135 

...1.4,649.... 

r-3. OBITER  2,PSIGMR/3 


0.0000.025  0.050  0.075  0.100  0.125  0.150  0.175  0.200  0.225  0.250  0.275  0.300 

time:  (SECONDS] 


Figure  C-14 

Timing  Discrepancies  Between  LX3965  Photo-Derived  (P},  and  Sensor-Derived  X 
Displacements  for  Head  (top)  and  T1  (bottom)  Anatomical  Origins 


Implementation  of  the  Maijoiization  Program 


F-3.0,ITE:R  2,PS!GI1fl  D/1,  P/50 


Legend 

LX390^  DRX50P  HOP  127  10.30g 

DMSQCb..  HOO 127  ...  1 0„ 

.  LX3901..  P.RXSOS . HOO 1 2.7 . 1.0..,.30g.. 


F-3. OBITER  2,PSrGMfi/3 


legend 

LX390^  DNX50P  H00I27  IO.3O9 
.LX390.i..DNX  S.OC  ....H0.0127 .  „ 
.kX.3901.JNX5p5 . HO.0.1 27 . 10- pOg ... 


TIME  (SECONDS) 


Figure  C-15 

Timing  Discrepancies  Between  LX3904  Photo-Derived  (P),  and  Sensor- 
Derived  X  Displacements  for  Head  (top)  and  T1  (bottom)  Anatomical  Origins 
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F- 

3.0, ITER  2,PSIGMfl  D/1,  P/100 

Legend 

LX3965  DflXSOP  HQQ135  H.64a 
LX3965  DRXSOC  H00135  14.6<lg 

LX3965  OflXSOS  H00135  14.64g 

-1.12S- 

“1.150- 

300 

to 

Cd 

^  “1.175- 

4 

QJ 

“1 .200- 

y 

CiJ 

-1.225- 

/ 

/ 

— 

LJ 

CJ 

01 

“1.250- 

// 

tn 

^  -1.275“ 

ij 

/ 

Q.OOO  0. 

025  0. 

F-3.( 

050  0 

3,  [IFF 

075  0. 

2, PS 

100  0. 

T 

IGMfi 

125  0. 

ME  (S 

3/1,  F 

ISO  Q 

ECOND 

/lOO, 

175  0. 

5) 

SLED 

200  0- 

nccEL 

225  0. 

-1 

250  0. 

275  0. 

Legend 

LX3965  DRX50P  H00135  14.B4q 

LX3965  DRXSOC  H00135  14.64g 

LX3965  DflXSOS  H00135  14.64g 

-1.100-] 

300 

-1.150- 

tn 

a: 

-1.175- 

4 

y 

CjJ 

z: 

4/ 

8— ' 

3: 

CJ 

2=  -1.225- 

- 

4^ 

O 

CE 

^  -1.250- 

h 

to 

^  -1.275- 

h 

- 1 . 300  - 

! 

- 1  - 

O.OOQ  0. 

D25  0. 

350  0. 

075  0.' 

00  0, 

TI 

25  0.' 

ME  (5f 

so  0. 
:C0ND5 

75  0.' 

1 

200  0. 

225  □. 

250  0.' 

275  0.' 

Figure  C-16 

Comparison  of  LX3965  Head  X  Displacements  for  Uncorrected  Marjorized  (top) 
and  Corrected  Marjorized  (bottom)  Reflecting  1  msec  Sled  Acceleration  Shift 
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P-3. OBITER  2,PSIBMR/3 


Legend 

LX3965  ONXSOP  H00135 

M.64q 

LX3965  DNXSOC 

H00135 

..Hr.eig.. 

LX3965  DNXSOS 

H00135 

M.64q 

_ U 

TIME  (SECONDS) 


F-3.0,ITER  2,PSIGMR/3,SLED  ACCEL  -1 


Legend 

LX3965 

DNXSOP 

HDD  135 

LX3965 

DNXSOC 

H00135 

.LX3965.. 

..QNxsgs... 

H00135 

TIME  (SECONDS) 


Figure  C-17 

Comparison  of  LX3965  Uncorrecte(d  Marjorized  (top)  and  Corrected  Marjorized 
(bottom)  T1  X  Displacements  Reflecting  1  msec  Sled  Acceleration  Shift 
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r 

-3.0, ITER  2,PSIGMR  D/1,  P/50 

Legand 

LX3904  DflXSOP  HQni27  ]0.30q 
LX3904  DRXSOCb  HOG  127  ]0.30q 

LX3904  DRXSOS  H00127  ]0.30g 

- 

300 

cr> 

cr. 

UJ 

s: 

-1 . 175 

/ 

X, 

E— • 

zz. 

LJ 

// 

LJ 

O 

CE 

-1.225- 

/ 

/ 

■•■X 

tn 

// 

/  / 

'' 

O.DQO  0, 

025  0. 

r-3. 

050  0 

OJTE 

Q75  0. 

R  2, PS 

100  0, 
T 

)IGMR 

- 1 - 

125  0.150  0. 

I  ME  (SECOND 

D/1,  P/50, 

175  0. 

5) 

SLED 

200  0- 

RCCEL 

225  0. 

-1 

250  0. 

275  0. 

Legend 

LX3904  DflXSOP  H00127  10.30a 

LX3904  DRXSOCb  HQ0127  10.30a 

.0X3904  .DRXSOS . HOQ 1 2.7 . 1.0,,„30g,.,. 

300 

cn 

cc 

ft!  -I.ISQ- 

LJ 

31 

f— ' 
z 

LJ 

S  -1.200- 

o 

cr 

—1  -!.22S- 

.4^ 

f 

"X 

cn 

^  -1.25G- 

/ 

/) 

J.DQO  0. 

D25  0. 

TSD  0. 

075  0. 

00  0. 
T1 

25  0. 

ME  (S 

so  0.' 
:C0ND5 

75  Q.' 

) 

JOO  0. 

>25  0.' 

250  0.' 

275  0. 

Figure  C-1  8 

Comparison  of  LX3904  Uncorrected  Marjorized  (top)  and  Corrected  Marjorized 
(bottom)  Head  X  Displacements  Reflecting  1  msec  Sled  Acceleration  Shift 
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r-3.0,IIER  2,PSIGMfl/3 


Legend 

LX390-!  DNX50P  HOP  1 27  1Q.30q 

.LX3901.DNX10C.„H001^^^^^ 
.L>:39,P.1.JNX505_,H0012^^^ . 1.0-3pg 


TIME  (SECONDS) 


P-3.0,  ITER  2,PSlGMfl/3,SLED  RCCEL  -1 


LX390-5 

Legend 

DNXSOP  H00127 

10.30g 

LX390^ 

DNXSOC 

H00127 

..J.Pr.lQg.. 

.LX390i. 

DNX50S 

H00127 

10.30g 

TIME  (SECONDS) 


Figure  C-19 

Comparison  of  LX3904  UncorrectecJ  Marjorized  (top)  and  Corrected  Marjorized 
(bottom)  T1  X  Displacements  Reflecting  1  msec  Sled  Acceleration  Shift 
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Appendix  D 
Tables 


Table  D-1 .  Data  Files  Used  by  the  Marjorization  Program 

Sample  User-Created  Data  Files 

Sample 
User-Created 
Data  File 
Name 

Input 

to 

Module 

Description 
of  Data 

Sample 

Data 

File 

Description 

of 

Data  File 

Name  Composition 

dllOl 

Module  1 

Table  D-5 

Appendix  B4 

1101  =Anat.Mt.iD 

d2201 

Module  1 

Table  D-5 

Appendix  B4 

2201  =Anat.Mt.lD 

data1 

Module  1 

Table  D-6 

Appendix  B5 

_ 

corrzero 

Module  1 

Table  D-11 

Appendix  B1 1 

filein 

Module  2 

Table  D-7 

Appendix  B6 

data2 

Module  2 

Table  D-8 

Appendix  B6 

_ 

oldc3959 

Module  2 

Table  D-11 

Appendix  BIO 

3959=  Sled  Run  No. 

dataS 

Module  3 

Table  D-9 

Appendix  B7 

— 

marjphodat 

Photo-Reformatting 

Table  D-10 

Appendix  B8 

_ 

pitmarj 

Plotting 

Table  D-1 1 

- 

— 

Sample  Program-Created  Data  Files 

Sample 

Program-Created 
Data  Files 

Output 

by 

Module 

Input 

to 

Sample 

Data 

File 

Description 

of 

Data 

head3959 

Photo-Reformatting  Module 

Module  1 

Appendix  B9 

Table  D-1 2 

neck3959 

Photo-Reformatting  Module 

Module  1 

_ 

Table  D-1 2 

out3959 

Module  1 

Module  2 

_ 

Page  9 

corr3959 

Module  2 

Module  3 

Appendix  BIO 

Table  D-11 
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Table  D-2.  E 

P 

valuation  of  Marjorized  Variables:’ 
hoto-Derived  Displacements. 

Standard  Factor  Set  Applied  to  Sigmas  for 

Head  Standard  Factor  Set:  1D,10P*** 

Sled 

Runs 

Head  Displacement  Variables 
(Variable  Definitions  in  Table  D-1 1) 

MAGNITUDE  OF 
CORRECTIONS** 

DAXSOC 

DAYSOC 

DAZSOC 

PHAOXC 

PHB02C 

PHC03C 

LX3852 

E 

E 

E 

E 

E 

E 

E 

LX3883 

E 

E 

E 

E 

E 

E 

E 

LX3893 

E 

E 

E 

E 

E 

E 

E 

LX3904 

G  + 

E 

G 

E 

E 

E 

E 

LX3924 

E 

E 

E 

E 

E 

E 

E 

LX3949 

E 

E 

E 

E 

E 

E 

E 

LX3957 

G 

E 

E 

E 

E 

E 

E 

LX3959 

E 

E 

E- 

E 

E 

E 

G 

LX3965 

G 

E 

E 

E 

E 

E- 

E 

LX3970 

G- 

E 

P 

E 

E 

E 

G 

LX3983 

E- 

E 

E 

E 

E 

E 

E 

LX3990 

E 

E 

E 

E 

E 

E 

E 

LX4552 

E 

E 

E 

E 

E 

E 

E 

Neck  Standard  Factor  Set:  3D,3P*** 

Sled  Runs 

Neck  Displacement  Variables 
(Variable  Definitions  in  Table  D-1 1 ) 

Magnitude  of 
Corrections** 

DNXSOC 

DNYSOC 

DNZSOC 

PNAOXC 

PNB02C 

PNC03C 

LX3852 

E- 

E 

E 

E 

E 

E 

G 

LX3883 

G 

E 

E 

E 

E 

E 

E 

LX3893 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

LX3904 

G  + 

E 

E 

E 

E 

E 

E 

LX3924 

E 

E- 

E 

E 

E 

E 

E 

LX3949 

G  + 

E 

E 

E 

E 

E 

E 

LX3957 

G 

E 

E 

E 

E 

E 

E 

LX3959 

G 

G  + 

E 

E 

E 

E 

E 

LX3965 

G 

E 

E 

E 

E 

E 

E 

LX3970 

G  + 

E 

E 

E 

E 

E 

E 

LX3983 

G 

E 

E 

E 

E 

E 

E- 

LX3990 

G  + 

E 

E 

E 

E 

E 

E 

LX4552 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 
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EVALUATION  CRITERIA: 


*Key  for  evaluating  Marjorized  variables  (differences  between  Marjorized  and  photo-derived  variables  —  peak  to  250 
msec). 

Description 

Key 

Linear 

Displacements 

Angular 

Displacements 

Excellent 

E 

=  1 ,25  cm 

<  =3  degrees 

Good 

G 

>1.25cm,  <=2.0cm 

>3  degrees,  <  =5  degrees 

Poor 

B 

>2. 0cm 

>  5  degrees 

**Key  for  evaluating  corrections  to  accelerometer  parameters  (overall  rating  is  that  of  the  component  with  the  largest 
error). 

Accelerometer 

Accelerometer 

Initial  Euler 

Description 

Key 

Orientation 

Sensitivity 

Angles 

Excellent 

E 

<  =3.5  degrees 

<  =3.5% 

<  =  3  degrees 

Good 

>3.5  degrees,  <  =5.5  degrees 

>3.5%,<=5.5% 

>3  degrees,  <  =  5  degrees 

Poor 

B 

>5.5  degrees 

>5.5% 

>5  degrees 

***Factor  Sets:  Factor  sets  are  defined  as  factors  dividing  variances  for  photo-derived  displacements. 
The  following  codes  were  used  to  identify  the  factor  sets: 

D- Displacement 
P  =  Euler  angles 
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Table  D-3.  Perturbation  Values  Used  for  Accelerometer  Parameters  in  Sample  Data  Files  d1 101 
and  (12201“ 

Line# 

Accelerometer 

Component  Perturbed 

Nominal 

Value 

2 

1 

Direction  Cosine  1 

0.0087178 

3 

1 

Direction  Cosine  2 

0.0087178 

4 

1 

Sensitivity 

0.0060000 

5 

2 

Direction  Cosine  1 

0.0061644 

6 

2 

Direction  Cosine  3 

0.0061644 

7 

2 

Sensitivity 

0.0060000 

8 

3 

Direction  Cosine  1 

0.0061644 

9 

3 

Direction  Cosine  3 

0.0061644 

10 

3 

Sensitivity 

0.0060000 

11 

4 

Direction  Cosine  1 

0.0061644 

12 

4 

Direction  Cosine  3 

0.0061644 

13 

4 

Sensitivity 

0.0060000 

14 

5 

Direction  Cosine  1 

0.0087178 

15 

5 

Direction  Cosine  2 

0.0087178 

16 

5 

Sensitivity 

0.0060000 

17 

6 

Direction  Cosine  1 

0.0087178 

18 

6 

Direction  Cosine  2 

0.0087178 

19 

6 

Sensitivity 

0.0060000 

20 

7 

Direction  Cosine  1 

0.0061644 

21 

7 

Direction  Cosine  3 

0.0061644 

22 

7 

Sensitivity 

0.0060000 

23 

8 

Direction  Cosine  1 

0.0061644 

24 

8 

Direction  Cosine  3 

0.0061644 

25 

8 

Sensitivity 

0.0060000 

26 

9 

Direction  Cosine  1 

0.0061644 

27 

9 

Direction  Cosine  3 

0.0061644 

28 

9 

Sensitivity 

0.0060000 

29 

Initial  Euler  Angle  1 

0.0435890 

30 

Initial  Euler  Angle  2 

0.0435890 

31 

Initial  Euler  Angle  3 

0.0435890 

"Perturbation  values  shown  in  Appendix  B4  for  the  data  files  dllOl  and  62201  for  anatomical  mounts  1101  and  2201 
respectively  are  3-cr  errors.  Errors  in  accelerometer  direction  cosines  and  sensitivities  were  based  on  multiple  calibrations  of  the 
T'Plate  and  were  suggested  by  an  earlier  study.^  Note  that  the  values  used  for  errors  in  the  orientation  of  the  accelerometers 
are  different  for  the  in-plane  accelerometers  and  those  perpendicular  to  the  plane  of  the  T-plate.  For  both  mounts  1 101  and 


2201,  accelerometers  1,  5,  and  6  were  perpendicular  to  the  plane  of  the  T-plate,  while  all  of  the  other  accelerometers  were 
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inplane.  3*a  errors  in  initial  Euler  angles  are  estimates  based  on  previous  studies  of  errors  in  photographic  data. 

Entries  in  d2201  for  anatomical  mount  2201  are  identical,  except  that  there  are  only  6  accelerometers  on  the  mount,  so  the  data 
file  only  contains  entries  for  accelerometers  1-6. 

Description  of  entries  in  data  file  dllOl  for  anatomical  mount  1101: 

Linel:  remarks  FORMAT  {A72) 

Lines  2-31:  perturbation  values  FORMAT{F1  5.0) 


Table  D-4.  Input  Variables  in  Data  File  dataV 

Line  1 : 

format(2i2) 

icorr 

Index  to  determine  if  corrections  are  to  be  applied  to  the  sensor  program  parameters 
to  generate  perturbed  values  for  the  follow 

Enter  0  if  corrections  are  not  to  be  applied. 

Enter  1  if  perturbations  are  to  be  applied. 

iprint 

Print  Index 

Enter  0  for  normal  output. 

Line  2: 

format(a6) 

test 

Sled  run  number 

Enter  6  character  sled  number,  e.g.,  LX3959. 

Line  3: 

format(a6) 

EOFEOF 

“Input  variables  for  Marjorization  Module  1  are  read  from  data  file  data1  (see  Appendix  B5) 
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Table  D-5.  Input  Variables  in  Data  File  filein' 

Line  1 : 

format(a1 1) 

filein 

The  name  of  the  file  which  contains  the  Marjorization  parameters  for  the  run. 
Enter  1 1  character  data  file  name,  e.g.,  data2 

Line  2: 

format(6f10.2) 

wth(i),i  =  1,6 

Factors  to  be  applied  to  the  denominators  of  the  a  priori  error  values  for  the 
head  photo-derived  displacements 

wth(1 )  =  factor  for  head  X  displacement,  nominal  value  1 .0 
wth(2)  =  factor  for  head  Y  displacement,  nominal  value  1 .0 
wth(3)  =  factor  for  head  Z  displacement,  nominal  value  1 .0 
wth(4)  =  factor  for  head  Euler  angle  A,  nominal  value  10.0 
wth(5)  =  factor  for  head  Euler  angle  B,  nominal  value  10.0 
wth(6)  =  factor  for  head  Euler  angle  C,  nominal  value  10.0 

Line  3: 

format(6f10.2) 

wtn(i),i  =  1,6 

Factors  to  be  applied  to  the  denominators  of  the  a  priori  error  values  for  the 
neck  photo-derived  displacements 

wtn(1)  =  factor  for  neck  X  displacement,  nominal  value  3.0 
wtn(2)  =  factor  for  neck  Y  displacement,  nominal  value  3.0 
wtn(3)  =  factor  for  neck  Z  displacement,  nominal  value  3.0 
wtn(4)  =  factor  for  neck  Euler  angle  A,  nominal  value  3.0 
wtn(5)  =  factor  for  neck  Euler  angle  B,  nominal  value  3.0 
vvtn(6)  =  factor  for  neck  Euler  angle  C,  nominal  value  3.0 

"Input  variables  for  Marjorization  Module  2  are  read  from  data  file  filein  (see  Appendices  B3  and  B6). 
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Table  D-6.  Input  Variables  in  Data  File  data2° 

Line  1 : 

format{a6,i2,i2,2f7.3) 

runid 

Sled  run  id,  e.g.,  LX3959 

iprint 

Print  index 

Set  to  0  for  normal  output 

iprior 

Index  to  determine  whether  a  priori  terms  will  be  used  to  constrain  the  solution  for 
the  corrections  to  the  sensor  program  parameters 

Set  to  0  for  unconstrained  solution. 

factor 

Factor  to  be  applied  to  a  priori  variances  for  accelerometer  direction  cosines  and 
sensitivities. 

Nominal  value  of  1  (1-10  normally  acceptable). 

factoreu 

Factor  to  be  applied  to  a  priori  variances  for  the  initial  Euler  angles. 

Nominal  value  of  1  (1-10  normally  acceptable). 

Line  2: 

format(7f7.4) 

seel  (i),i  =  1,7 

Endpoint  of  the  ith  time  interval  during  the  sled  run,  where  the  time  window  of 
interest  is  subdivided  into  7  intervals;  used  for  head  calculations. 

Nominal  values  shown  in  Appendix  B3. 

Line  3: 

format(7f7.4)  4 

sec2(i),i  =  1,7 

Endpoint  of  the  ith  time  interval  during  the  sled  run,  where  the  time  window  of 
interest  is  subdivided  into  7  intervals;  used  for  neck  calculations 

Nominal  values  shown  in  Appendix  B3. 

Line  4: 

format(7i3) 

iskp1{i),i  =  1,7 

Determines  the  number  of  data  points  used  in  calculations  from  the  ith  time 
interval.  The  value  selected  for  iskpl  (i)  causes  iskpl  (i)-1  points  to  be  skipped  in 
the  ith  interval;  used  for  head  calculations. 

Nominal  values  shown  in  Appendix  B3. 

Line  5: 

format(7i3) 

iskp2(i),i  =  1,7 

Determines  the  number  of  data  points  used  in  calculations  from  the  ith  time 
interval.  The  value  selected  for  iskp2(i)  causes  iskp2(i)-1  points  to  be  skipped  in 
the  ith  interval;  used  for  neck  calculations. 

Nominal  values  shown  in  Appendix  B3. 

•Input  variables  for  Marjorization  Module  2  are  read  from  data  file  data2  (see  Appendices  B3  and  B6) 
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Table  D-7.  Input  Variables  in  Data  File  data3‘ 

Line  1: 

format(2i2) 

icorr 

Index  to  determine  if  corrections  are  to  be  applied  to  the  sensor  program  parameters  to 
obtain  Marjorized  variables. 

Enter  0  if  corrections  are  not  to  be  applied. 

Enter  1  if  Marjorized  variables  are  to  be  obtained  by  applying  corrections  to  the  sensor 
program  parameters. 

iprint 

Print  Index 

Enter  0  for  normal  output. 

Line  2: 

format(a6) 

test 

Sled  run  number 

Enter  6  character  sled  number,  e.g.,  LX3959. 

Line  3: 

format(a6) 

EOFEOF 

■Input  variables  for  Marjorization  Module  3  are  read  from  data  file  data3  (see  Appendix  B7) 


Table  D-8.  Input  Variables  in  Data  File  marjphodaf 

Line  1 : 

format(i2,1x,a6) 

iunit 

The  unit  (file)  containing  the  EZFLOW  reduced  photo  variables. 

The  nominal  value  is  7,  since  this  is  the  unit  on  which  EZFLOW  expects  the  reduced 
photo  variables  to  be  stored. 

runid 

Sled  run  id,  e.g.,  LX3959 

■Input  variables  for  the  Marjorization  photo-reformatting  module,  marjphopgm.f,  are  read  from  data  file  marjphodat  (see  Appendix 
B8). 
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Table  D-9. 

Variables  Included  in  Plot  Specification  Module  pitmarj 

PLOT  1 : 

DAXSOP,  DAXSOS,  DAXSOC 

The  component  of  linear  displacement  of  the  head  anatomical  origin  along  the  X  axis 
of  the  sled  coordinate  system  with  respect  to  the  sled  coordinate  system.  The  sled 
coordinate  system  is  aligned  with  the  laboratory  coordinate  system  and  translates 
with  the  sled. 

PLOT  2: 

DAYSOP,  DAYSOS,  DAYSOC 

The  component  of  linear  displacement  of  the  head  anatomical  origin  along  the  Y  axis 
of  the  sled  coordinate  system  with  respect  to  the  sled  coordinate  system.  The  sled 
coordinate  system  is  aligned  with  the  laboratory  coordinate  system  and  translates 
with  the  sled. 

PLOT  3; 

DAZSOP,  DAZSOS,  DAZSOC 

The  component  of  linear  displacement  of  the  head  anatomical  origin  along  the  Z  axis 
of  the  sled  coordinate  system  with  respect  to  the  sled  coordinate  system.  The  sled 
coordinate  system  is  aligned  with  the  laboratory  coordinate  system  and  translates 
with  the  sled. 

PLOT  4: 

PHAOXP,  PHAOXS,  PHAOXC 

Angle  of  rotation  of  the  head  about  the  X  axis  of  the  head  anatomical  coordinate 
system  as  derived  from  mouth-mount  data.  The  head  anatomical  coordinate  system 
is  initially  aligned  with  the  laboratory  coordinate  system 

PLOT  5: 

PHB02P,  PHB02S,  PHB02C 

Same  as  above  except  about  the  carried  Y  axis. 

PLOT  6: 

PHC03P,  PHC03S,  PHC03C 

Same  as  above  except  about  the  carried  Z  axis. 

PLOT  7: 

DNXSOP,  DNXSOS,  DNXSOC 

The  component  of  linear  displacement  of  the  T1  anatomical  origin  along  the  X  axis  of 
the  sled  coordinate  system  with  respect  to  the  sled  coordinate  system.  The  sled 
coordinate  system  is  aligned  with  the  laboratory  coordinate  system  and  translates 
with  the  sled. 

PLOT  8: 

DNYSOP,  DNYSOS,  DNYSOC 

The  component  of  linear  displacement  of  the  T1  anatomical  origin  along  the  Y  axis  of 
the  sled  coordinate  system  with  respect  to  the  sled  coordinate  system.  The  sled 
coordinate  system  is  aligned  with  the  laboratory  coordinate  system  and  translates 
with  the  sled. 

PLOT  9: 

( 

DNZSOP,  DNZSOS,  DNZSOC 

The  component  of  linear  displacement  of  the  T1  anatomical  origin  along  the  Z  axis  of 
the  sled  coordinate  system  with  respect  to  the  sled  coordinate  system.  The  sled 
coordinate  system  is  aligned  with  the  laboratory  coordinate  system  and  translates 
with  the  sled. 

PLOT  10: 

> 

i 

PNAOXP,  PNAOXS,  PNAOXC 

^ngle  of  rotation  of  the  head  about  the  X  axis  of  the  T1  anatomical  coordinate 
system  as  derived  from  T1 -mount  data.  The  T1  anatomical  coordinate  system  is 
nitially  alignod  with  the  laboratory  coordinate  system. 
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Table  D-9.  Variables  Included  in  Plot  Specification  Module  pitman 


PLOT  1 1 :  PNB02P,  PNB02S,  PNB02C 
_ Same  as  above  except  about  the  carried  Y  axis. 


PLOT  12:  PNC03P,  PNC03S,  PNC03C 

Same  as  above  except  about  the  carried  Z  axis. 


Table  D-10.  Description  of  Data  in  Reformatted  Photo  Data  Files" 


_  Line  1 :  format(a4,8(10x,a6)) 


head  variable  names"  "time,"  "daxsop,"  "daysop,"  "dazsop,"  "4hoo1p,"  "4hoo2p,"  "4hoo3p," 
_ _ "4hoo4p" 


neck  variable  names"  "time,"  "dnxsop,"  "dnysop,"  "dnzsop,"  "4noo1p,"  "4noo2p,"  "4noo3p," 
_ "4noo4p" _ 


_  Lines  2  ...  End  of  File:  format(e14.8,7(1x,e15.8)) 


Each  line  contains  values  for  the  8  variables  input  on  Line  1 . _ 

•The  Marjorization  Photo  Reformatting  Module  reformats  the  photo-derived  linear  displacement  variables  and  quaternions  in  the 
format  required  by  Marjorization  module  2.  These  data  are  output  in  the  following  format  (see  Appendix  B9): 

"See  Table  D-9  for  a  definition  of  the  variables  beginning  with  "d",  the  photo  linear  displacement  variables.  The  variables 
beginning  with  "4h"  are  head  quaternions  and  those  beginning  with  "4n"  are  neck  quaternions,  where  quaternions  are  the  four 
variables  which  define  the  angular  orientation  of  the  head  or  T1  anatomical  coordinate  system  relative  to  the  laboratory 
coordinate  system. 
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Table  D-1 1 .  Description  of  the  Output  Solution  File  from  Marjorization  Module  2  and  Other 
_ Correction  Files* _ 

The  solution  output  by  Marjorization  Module  2  is  written  to  a  file  whose  name  is  "corr"  concatenated  with  the  four-digit 
numeric  portion  of  the  sled  run  number,  e.g.,  corr3959  (Appendix  BIO).  This  file  is  read  by  Marjorization  Module  1  for 
successive  iterations  and  must  be  initialized  to  all  zero  components  prior  to  the  execution  of  Module  1  for  the  first 
iteration.  This  initialization  is  accomplished  in  the  shell  script  xqtmarj  (Appendix  B1)  by  providing  the  file  corrzero 
(Appendix  B1 1)  with  all  zero  error  components,  which  is  copied  to  the  file  corr3959  prior  to  the  first  iteration.  The  file 
o!dc3959  containing  corrections  from  the  previous  iteration  is  also  required  as  input  to  Marjorization  Module  2,  and 
corr3959  must  not  exist  when  the  execution  of  Module  2  is  initiated.  Script  run2  (Appendix  B2)  takes  care  of  both  of 
these  requirements  automatically  by  copying  corr3959  into  oldc3959  and  deleting  corr3959  prior  to  executing 
Marjorization  module  2.  The  file  corr3959  generated  by  Module  2  for  the  final  iteration  is  input  to  Marjorization  Module  3, 
which  generates  the  Marjorized  variables. 


The  error  components  contained  in  corr3959  are  presented  in  a  more  readable  form  in  the  output  from  Module  3  (Table 
D-1 4),  in  which  the  error  components  are  shown  as  the  orientation  error  for  each  accelerometer  in  degrees  and  the 
sensitivity  error  as  a  percentage.  Initial  Euler  angle  errors  are  also  presented  in  degrees. 


Line  1 : 

Identifying  information  for  the  first  anatomical  mount:  Run  No.,  data  type,  and 
number  of  correction  components.  In  the  example  shown  in  Appendix  BIO,  the 
first  mount  has  9  accelerometers  and  therefore  30  error  components. 

Line  2: 

Correction  components  1-8,  where  the  computed  correction  components  are  for 
the  indicated  accelerometer  error  component  for  the  first  mount: 

Correction  Component 

Accelerometer 

Component 

1 

1 

Direction  Cosine  1 

2 

1 

Direction  Cosine  2 

3 

1 

Sensitivity 

4 

2 

Direction  Cosine  1 

5 

2 

Direction  Cosine  3 

6 

2 

Sensitivity 

7 

3 

Direction  Cosine  1 

8 

3 

Direction  Cosine  3 

Line  3: 

Correction  components  9-16,  where  the  computed  correction  components  are  for 
the  indicated  accelerometer  error  component  for  the  first  mount: 

Correction  Component 

Accelerometer 

Component 

9 

3 

Sensitivity 

10 

4 

Direction  Cosine  1 

11 

4 

Direction  Cosine  3 

12 

4 

Sensitivity 

13 

5 

Direction  Cosine  1 

14 

5 

Direction  Cosine  2 

15 

5 

Sensitivity 

16 

6 

Direction  Cosine  1 
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Table  D-1 1 .  Description  of  the  Output  Solution  File  from  Marjorization  Module  2  and  Other 
Correction  Files® 


Correction  components  1 7-24,  where  the  computed  correction  components  are 
for  the  indicated  accelerometer  error  component  for  the  first  mount: 


Correction  Component 


17 


Accelerometer 


Component 


Direction  Cosine  2 


Sensitivity 


Direction  Cosine  1 


Direction  Cosine  3 


Sensitivity 


Direction  Cosine  1 


Direction  Cosine  3 


Sensitivity 


Correction  components  25-30,  where  the  computed  correction  components  are 
for  the  indicated  accelerometer  error  component  for  the  first  mount: 


Correction  Component 


25 


Accelerometer 


Component 


Direction  Cosine  1 


Direction  Cosine  3 


Sensitivity 


Initial  Euler  Angle  1 


Initial  Euler  Angle  2 


Initial  Euler  Angle  3 


Identifying  Information  for  the  second  anatomical  mount:  Run  No.,  data  type,  and 
number  of  correction  components.  In  the  example  shown  in  Appendix  BIO  the 
second  mount  has  6  accelerometers  and  therefore  21  error  components. 


Correction  components  1-8,  where  the  computed  correction  components  are  for 
the  accelerometer  error  component  for  the  second  mount  described  in  Line  2 
above. _ 


Correction  components  9-16,  where  the  computed  correction  components  are  for 
the  accelerometer  error  component  for  the  second  mount  described  in  Line  3 
above. _ 


Correction  components  17-21,  where  the  computed  correction  components  are 
for  the  Indicated  accelerometer  error  component  for  the  second  mount: 


Correction  Component 


Accelerometer 


Component 


Direction  Cosine  2 


Sensitivity 


Initial  Euler  Angle  1 


Initial  Euler  Angle  2 


Initial  Euler  Angle  3 


•The  following  variable  descriptions  apply  to  the  solution  output  by  Marjorization  module  2  as  shown  in  the  sample  output  file 
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corr3959  in  Appendix  BIO.  These  descriptions  also  apply  to  files  corrzero  and  oldc3959. 


Table  D-12.  Evaluation  of  Marjorized  Variables® 


Customized  Factor  Sets  Applied  to  Sigmas  for  Photo-Derived  Displacements 


Sled 

Runs 

Factor 

Set" 

Head  displacement  variables  (definitions  in  Table  D-1 1) 

Magnitude  of 
Corrections'" 

DAXSOC 

DAYSOC 

DAZSOC 

PHAOXC 

PHB02C 

PHC03C 

LX3904 

3D,50P 

G  + 

E 

E 

E 

E 

E 

G  + 

LX3970 

3D&P,10Z 

E 

E 

E 

E 

E 

E 

G 

Sled 

Runs 

Factor 

Set" 

Neck  displacement  variables  (definitions  in  Table  D-1 1) 

Magnitude  of 
Corrections'" 

DNXSOC 

DNYSOC 

DNZSOC 

PNAOXC 

PNB02C 

PNC03C 

LX3959 

10D,10P 

G  + 

E 

E 

E 

E 

E 

E 

The  following  codes  were  used  to  identify  the  factor  sets: 
D  =  Displacement 
P  =  Euler  angles 
Z  =  Z  displacement 

Evaluation  criteria: 


“Key  for  evaluating  marjorized  variables 

(Differences  between  marjorized  &  photo-derived  variables  —  peak  to  250  msec) 


Description 

Key 

Linear 

Displacements 

Angular 

Displacements 

Excellent 

E 

<  =  1 .25  cm 

<  =3  degrees 

Good 

G 

>  1 .25  cm,  <  —  2.0  cm 

>3  degrees,  <  =  5  degrees 

Poor 

P 

>2.0  cm 

>5  degrees 

**Factor  sets: 

Factor  sets  are  defined  as  factors  dividing  variances  for  photo-derived  displacements. 


‘'Key  for  evalual 
(Overall  rating  i 

ing  CO 

s  that 

rrections  to  accelerometer  parameters 
)f  the  component  with  the  largest  error) 

Description 

Key 

Accelerometer 

Orientation 

Accelerometer 

Sensitivity 

Initial  Euler 

Angles 

Excellent 

E 

<  -3.5  degrees 

<  =3.5% 

<  =  3  degrees 

Good 

>3.5  degrees,  <  =5.5  degrees 

>3.5%,<  =5.5% 

>3  degrees,  <  =5  degrees 

Poor 

B 

>5.5  degrees 

>5.5% 

>5  degrees 
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Table  D-13.  Evaluation  of  Marjorized  Variables® 


Sled 

Runs 

Head  displacement  variables  (definitions  in  Table  D-11) 

Magnitude  of 
Corrections' 

I  acior 

Set" 

DAXSOC 

DAYSOC 

DAZSOC 

PHAOXC 

PHB02C 

PHC03C 

LX3904 

1D,20P 

E- 

E 

E 

E 

E 

E 

E- 

LX3965 

ID,  loop 

E 

E 

E 

E 

E 

G  + 

E 

Sled 

Runs 

Fector 

Neck  displacement  variables  (definitions  in  Table  D-11) 

Magnitude  of 
Corrections*^ 

Set" 

DNXSOC 

DNYSOC 

DNZSOC 

PNAOXC 

PNB02C 

PNC03C 

LX3904 

3D,3P 

E 

E 

E 

E 

E 

E 

E 

LX3965 

3D,3P 

E 

E 

E 

E 

E 

E 

E 

Customized  factor  sets  applied  to  sigmas  for  photo-derived  displacements  sled  acceleration;  trace  shifted  to  the  left  by  1  msec 
The  following  codes  were  used  to  identify  the  factor  sets: 

D  =  Displacement 
P  =  Euler  angles 
H  =  Head 
N^Neck 

Z  =  Z  displacement 
EVALUATION  CRITERIA: 


"Key  for  evaluating  Marjorized  variables 

(Differences  between  Marjorized  &  photo-derived  variables  —  peak  to  250  msec) 


Description 

Key 

Linear 

Displacements 

Angular 

Displacements 

Excellent 

E 

<  =  1 .25  cm 

<  ~3  degrees 

Good 

G 

>  1 .25  cm,  <  =2.0  cm 

>3  degrees,  <  =  5  degrees 

Poor 

P 

>2.0  cm 

>5  degrees 

‘Tactor  sets: 

Factor  sets  are  defined  as  factors  dividing  variances  for  photo-derived  displacements. 

Note  that  the  old  weighting  scheme  using  up  to  60  data  points  was  used  for  these  runs,  since  this  work  was  done  before  the 
new  weighting  scheme  which  uses  up  to  150  data  points  was  implemented. 
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’'Key  for  evaluat 
(Overall  rating  is 

ng  cor 

that  c 

rections  to  accelerometer  parameters 
f  the  component  with  the  largest  error) 

Description 

Key 

Accelerometer 

Orientation 

Accelerometer 

Sensitivity 

Initial  Euler 

Angles 

Excellent 

E 

<  =3.5  degrees 

<  =3.5% 

<  =  3  degrees 

Good 

G 

>3.5  degrees,  <  =5.5  degrees 

>3.5%,  <  =5.5% 

>3  degrees,  <  =5  degrees 

Poor 

P 

>  5.5  degrees 

>5.5% 

>5  degrees 

Table  D-14.  Computed  Corrections  on  Second  Iteration  —  LX3959” 

Head  corrections:  baseline  factors  applied  to  denominators  of  head  photo-derived  displacement 
variances  (linear  displacement  factor:  1,  angular  displacement  factor:  10) 

Accelerometer 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Orientation  (degrees) 

1.81 

1.11 

0.46 

1.94 

3.36 

2.24 

0.55 

0.79 

0.90 

Sensitivity  (%) 

-3.28 

1.37 

0.04 

0.81 

-3.50 

-3.25 

0.54 

0.40 

0.16 

Euler  angle 

1 

2 

3 

Initial  orientation 
(degrees): 

0.63 

-0.98 

-0.48 

Neck  corrections:  baseline  factors  applied  to  denominators  of  neck  photo-derived  dis| 
variances  (linear  displacement  factor:  5,  angular  displacement  factor:  5) 

Diaceme 

mt 

Accelerometer 

1 

2 

3 

4 

5 

6 

Orientation  (degrees) 

1.60 

4.75 

0.90 

1.70 

1.03 

1.61 

Sensitivity  (%) 

2.69 

-2.10 

0.26 

2.25 

-0.02 

-2.34 

Euler  angle 

1 

2 

3 

Initial  orientation 
(degrees): 

0.38 

-1.06 

-0.67 

'Note  that  the  old  weighting  scheme  using  up  to  60  data  points  was  used  for  these  runs,  since  this 
work  was  done  before  the  new  weighting  scheme  which  uses  up  to  1  50  data  points  was  implemented. 
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